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1.  INTRODUCTION: 


Prostate  cancer  remains  the  most  common  non-cutaneous  malignancy  in  men  in  USA.  Active  surveillance 
can  be  applied  to  most  patients.  However,  around  15%  of  prostate  cancer  patients  are  diagnosed  with  high 
risk  disease,  and  these  patients  have  significantly  higher  chances  of  developing  systemic  metastasis  or  local 
recurrence,  therefore  more  intensive  treatment  is  required.  Treatment  options,  including  surgery,  radiation 
therapy,  and  androgen-deprivation  therapy,  or  their  combinations,  are  often  ineffective  in  a  subset  of  high 
risk  patients.  Most  patients  with  high-risk  disease  relapse  after  therapy,  and  there  are  no  effective  preventive 
measures.  Immunotherapy  holds  great  potential  in  treating  prostate  cancer.  However,  major  challenges 
remain.  We  propose  an  innovative  immunotherapy  approach  to  kill  tumor  cells,  one  that  combines  the 
specificity  of  gene  circuits  with  the  ability  of  Surface  T-cell  Engagers  (STEs).  STEs  can  trigger  antigen- 
independent  T-cell  killing  of  tumor  cells  in  vitro  and  in  vivo.  This  approach,  Tumor  Recruited 
Immunocellular  Therapy  (STRICT)  is  a  unique  and  innovative  combination  of  cancer-detecting  gene 
circuits  with  anti-cancer  immunotherapies.  Synthetic  gene  circuits,  administered  systemically,  will  be 
selectively  activated  only  in  cancer  cells,  resulting  in  the  surface  display  of  STE  and  the  secretion  of  potent 
immunomodulatory  molecules.  The  STEs  will  recruit  T  cells  to  target  and  kill  the  tumor.  Combining  the 
advantages  of  localized  immunotherapy  with  those  of  systemic  delivery,  STRICT  would  obviate  the  need 
for  systemic  chemotherapy  and  protect  against  tumor  relapse.  In  addition,  STRICT  is  expected  to  overcome 
the  key  obstacles  associated  with  other  therapies  currently  available  for  prostate  cancer. 

2.  KEYWORDS:  synthetic  gene  circuit;  cancer  immunotherapy,  prostate  cancer 

3.  ACCOMPLISHMENTS: 

What  were  the  major  goals  of  the  project? 


Months 

Percentage  of 
Completion 

Major  Task  1:  Create,  optimize,  and  validate 
cancer-detecting  circuits  that  display  Surface  T- 
cell  Engagers  (STEs)  from  within  tumors. 

1-9 

75% 

Major  Task  2:  Identify  the  minimal  percentage  of 
tumor  cells  that  need  to  be  targeted  by  TIGRIS 
for  in  vivo  efficacy. 

9-15 

0% 

Major  Task  3:  Determine  whether  TIGRIS  can 
eliminate  metastatic  tumors  via  systemic  delivery. 

3-36 

0% 

Major  Task  4:  Build  synthetic  safety  mechanisms 
for  TIGRIS  circuits. 

15-30 

0% 

Major  Task  5:  Test  if  TIGRIS  can  elicit  immune 
memory  to  protect  against  future  tumor  relapse. 

15-36 

0% 

What  was  accomplished  under  these  goals? 
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The  major  activities  throughout  this  report  period  focused  on  building  a  versatile  Tumor  Immunotherapy 
by  Gene-circuit  Recruited  Immunomodulatory  Systems  (TIGRIS)  for  prostate  cancer.  In  this  research 
period,  by  using  ovarian  cancer  model,  we  have  validated  that  TIGRIS  have  the  potential  to  trigger  robust 
immunotherapy  in  vitro  and  in  vivo  (see  appendix  for  an  accepted  manuscript).  We  have  optimized  the  key 
circuit  building  blocks  for  creating  highly  effective  RNA-based  synthetic  gene  circuit  for  tumor-specific 
immunomodulation.  We  also  optimized  circuit  architecture  that  enables  this  system  to  be  high  modular  for 
easy  adaptation  for  treating  other  cancer  type  and  for  future  clinical  translation.  In  addition,  we  have  created 
a  method  to  build  synthetic  promoters  that  perform  better  than  native  promoters  for  targeting  cancers 
(Figure  1).  These  promoters  directly  serve  as  sensors  for  TIGRIS  system  and  we  further  extend  this  method 
to  build  a  library  of  synthetic  promoters  (  >  6000  different  synthetic  promoters)  that  serve  as  a  sensor  library 
for  TIGRIS  system.  This  library  can  be  used  to  detect  multiple  different  cancer  types.  Currently,  we  are  able 
to  find  new  sensors  for  building  this  system  toward  treating  prostate  cancer  (Table  1).  The  next  step  forward 
is  to  validate  whether  combining  the  sensors  with  the  optimal  circuit  architecture,  the  system  remains  highly 


specific  to  prostate  cancer  cells  but  not  to  other  major  tissue  types  in  vitro. 
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Figure  1  -  Specificity  comparison  of  natural  and  synthetic  promoters.  (A)  Natural  occurring  promoters  may 
have  hundreds  of  possible  TFBS.  In  the  ovarian  specific  promoter  1  (OSPl)  there  are  more  than  a  hundred 
possible  TFBS,  each  color  represents  a  different  TFBS  family.  (B)  Synthetic  promoter  for  Jun  (API)  contains 
several  TFBS  for  Jun  upstream  to  minimal  adenovirus  promoter  (lateADEp)  and  TSS.  (C)  the  synthetic  Jun(APl) 
we  constructed  is  activated  specifically  in  one  cell  line  (HCT)  and  not  the  others.  OSPl  appears  nonspecific, 
being  activated  in  most  cell  lines  shown.  Both  promoters  show  comparable  activation  levels  in  multiple  cancer 
(HCT,  MDA-453,  MDA-231,  OVCAR8)  and  normal  cell  lines  (MCFIOA,  MCF12A,  adult  HDF). 


What  opportunities  for  training  and  professional  development  has  the  project  provided? 

This  project  provides  extensive  training  for  postdoe  fellow  Dr.  Ming-Ru  Wu  on  designing  eaneer- 
deteetion  eireuits,  optimizing  eireuit  inputs  and  outputs,  and  eaneer  immunotherapy.  It  also  provided  Dr. 
Ming-Ru  Wu  with  hands  on  experienee  on  working  with  prostate  eaneer  models. 

How  were  the  results  disseminated  to  communities  of  interest? 
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Nothing  to  report. 


What  do  you  plan  to  do  during  the  next  reporting  period  to  accomplish  the  goals? 

We  will  further  utilize  the  identified  sensors  to  build  TIGRIS  system  for  prostate  eaneer.  We  will  validate 
the  system  is  specific  to  prostate  cancer  cells  but  not  normal  tissues.  In  addition,  we  will  identify  the 
minimal  percentage  of  tumor  cells  that  need  to  be  targeted  by  TIGRIS  system  to  mediate  in  vivo  efficacy. 
Furthermore,  we  will  determine  whether  TIGRIS  can  eliminate  metastatic  tumors  via  systemic  circuit 
delivery. 


promoter 

Transcription  factor  binding  site 

MFI 

(PCS) 

Fold  to 
background 
MFI 

MFI 

(DU145) 

Fold  to 
background 
MFI 

1 

SPDEFmus 

775 

16 

513 

11 

2 

MYC 

46 

1 

22 

0 

3 

MYC 

49 

1 

44 

1 

4 

MYC/USFl 

47 

1 

44 

1 

5 

MY  CN_mus/CLK_d.melanogaster 

47 

1 

41 

1 

6 

TJ  d.melanogaster  (MAFB) 
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6 

131 

3 

7 

SPDEFmus 

511 

10 

588 

12 

8 

FOXA 

47 

1 

41 

1 

9 

MAFK 

148 

3 

80 

2 

10 

FOXA3 

416 

9 

86 

2 

11 

E2F1/E2F4 

46 

1 

43 

1 

12 

E2F 

52 

1 

83 

2 

13 

MAFG 

60 

1 

59 

1 

14 

MAFK 

167 

3 

86 

2 

Table  1:  Identification  of  promoters  that  are  active  in  prostate  cancers.  The  activity  of  14  breast 
cancer  specific  promoters  were  tested  on  2  prostate  cancer  cell  lines  (PCS,  DU145)  by  lentiviral 
transduction.  MFI  =  mean  fluorescence  intensity.  Fold  to  background  MFI  was  calculated  by  dividing 
the  MFI  of  individual  promoters  to  the  MFI  of  infection  control  cells.  2-3  promoters  exhibit  >  9  fold 
higher  activity  than  background  fluorescence  and  can  be  further  used  to  build  the  TIGRIS  system. 

4.  IMPACT: 

What  was  the  impact  on  the  development  of  the  principal  discipline(s)  of  the  project? 

Our  data  utilizing  TIGRIS  system  for  treating  ovarian  cancer  had  demonstrated  the  feasibility  and  the 
potential  robust  efficacy  of  TIGRIS  system.  If  the  adaptation  of  TIGRIS  system  to  prostate  cancer  is 
successful,  this  technology  would  represent  a  major  step  forward  in  our  ability  to  eliminate  prostate  cancer, 
especially  metastatic  and  recurrent  disease,  in  a  scalable  and  programmable  fashion. 

What  was  the  impact  on  other  disciplines? 

Our  data  highlighted  the  potential  power  of  utilizing  gene  circuit  to  treat  unmet  clinical  needs. 
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What  was  the  impact  on  technology  transfer? 

Nothing  to  report. 

What  was  the  impact  on  society  beyond  science  and  technology? 

Nothing  to  report. 

5.  CHANGES/PROBLEMS 

Changes  in  approach  and  reasons  for  change 

Nothing  to  report. 

Actual  or  anticipated  problems  or  delays  and  actions  or  plans  to  resolve  them 

Nothing  to  report. 

Changes  that  had  a  significant  impact  on  expenditures 

Nothing  to  report. 

Significant  changes  in  use  or  care  of  human  subjects,  vertebrate  animals,  biohazards,  and/or  select 
agents 

Nothing  to  report. 

Significant  changes  in  use  or  care  of  human  subjects 

Nothing  to  report. 

Significant  changes  in  use  or  care  of  vertebrate  animals 

Nothing  to  report. 

Significant  changes  in  use  of  biohazards  and/or  select  agents 
Nothing  to  report. 

6.  PRODUCTS: 

Publications,  conference  papers,  and  presentations 
Journal  publications. 

Nissim  L,  Wu  MR  (Co-first  author),  Pery  E,  Binder-Nissim  A,  Suzuki  HI,  Stupp  D,  Wehrspaun  C, 
Tabach  Y,  Sharp  PA,  and  Lu  TK.  2017.  Synthetic  RNA-based  immunomodulatory  gene  circuits  for 
cancer  immunotherapy.  Cell  (accepted,  awaiting  publication);  acknowledgement  of  federal  support  (yes). 

Books  or  other  non-periodical,  one-time  publications. 

Nothing  to  report. 
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Other  publications,  conference  papers,  and  presentations. 

Nothing  to  report. 

Website(s)  or  other  Internet  site(s) 

Nothing  to  report. 

Technologies  or  techniques 

Nothing  to  report. 

Inventions,  patent  applications,  and/or  licenses 

Patent  submitted  by  Massachusetts  Institute  of  Technology:  Lu  et  al.  (2016)  U.S.  TUMOR 
IMMUNOTHERAPY.  62/181 906,M0656.70356 


Other  Products 

Nothing  to  report. 


7.  PARTICIPANTS  &  OTHER  COLLABORATING  ORGANIZATIONS 
What  individuals  have  worked  on  the  project? 


Name: 

Ming-Ru  Wu 

Project  Role: 

Postdoc  Researcher 

Researcher  Identifier  (e.g. 
ORCID  ID): 

orcid.org/0000-0002-2533-8484 

Nearest  person  month 
worked: 

12 

Contribution  to  Project: 

Dr.  Wu  has  been  performing  designing  and  optimizing 
TIGRIS  system  for  prostate  cancer 

Funding  Support: 

Department  of  Defense  (W81XWH-16-1-0452). 

Has  there  been  a  change  in  the  active  other  support  of  the  PD/PI(s)  or  senior/key  personnel  since  the 
last  reporting  period? 

Nothing  to  report. 

What  other  organizations  were  involved  as  partners? 

Nothing  to  report. 

8.  SPECIAL  REPORTING  REQUIREMENTS 

Nothing  to  report. 
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9.  APPENDICES: 


Journal  publications. 

Nissim  L,  Wu  MR  (Co-first  author),  Pery  E,  Binder-Nissim  A,  Suzuki  HI,  Stupp  D,  Wehrspaun  C,  Tabach 
Y,  Sharp  PA,  and  Lu  TK.  2017.  Synthetic  RNA-based  immunomodulatory  gene  circuits  for  cancer 
immunotherapy.  Cell  (accepted,  awaiting  publication);  acknowledgement  of  federal  support  (yes). 
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Dear  Tim 


It  is  a  pleasure  to  let  you  know  that  we  will  be  publishing  your  paper  and  that  it  is  being  tentatively  scheduled  for  the 
November  30  issue  of  Cell.  Congratulations  on  a  very  interesting  story;  we  are  pleased  to  be  bringing  it  to  the  attention 
of  our  readers.  In  the  next  week  or  so,  our  Production  department  will  contact  you  with  your  online  publication  date  that 
will  be  ahead  of  the  print  issue  date. 

At  this  time,  I  also  wanted  to  introduce  our  Managing  Editor,  Trina  Arpin,  who  heads  up  the  production  team  that  will 
guide  your  paper  through  copyediting,  proofing,  and  publication.  They  will  be  contacting  you  shortly  with  some  initial 
information  on  the  production  process.  If  you  anticipate  being  unreachable  by  email  at  any  time  between  now  and  the 
publication  date,  please  send  them  (tarpin@cell.com)  the  email  address  of  an  alternate  author  who  can  rapidly  respond 
to  queries  from  our  production  department. 

In  the  coming  weeks,  your  assigned  copyeditor  will  be  in  contact  with  you  concerning  queries  on  the  manuscript  and 
PDF  proofs.  We  would  greatly  appreciate  it  if  you  can  keep  an  eye  out  for  emails  concerning  these  finishing  steps  so  that 
we  can  move  your  paper  efficiently  through  to  publication.  Please  note  that  at  this  stage  we  can  only  allow  minor 
corrective  changes  to  the  manuscript,  and  none  that  will  add  or  change  results  or  significantly  alter  the  length/style  of 
the  title,  abstract  or  manuscript  text.  If  any  concerns  about  this  arise,  please  contact  me  directly.  Please  see  below  for 
further  details  about  proofs. 

A  note  about  press,  our  general  policy  with  regard  to  press  coverage  is  that  authors  are  free  to  talk  with  the  popular 
press  starting  one  week  before  publication.  Your  press  office  may  issue  its  own  press  release  provided  that  it  is 
embargoed  until  12:00pm  NOON  EST  on  the  day  of  online  publication  (as  mentioned  above,  we'll  let  you  know  about  the 
online  publication  date  soon).  We  also  allow  authors  to  discuss  their  work  in  press  with  other  scientific  journals.  Because 
journals  could  in  principle  use  that  information  to  speed  publication  of  their  own  competitive  manuscripts,  it  is  best  if 
you  refer  any  inquiries  to  us.  Once  the  requesting  journal  has  confirmed  that  they  do  not  have  any  related  work  under 
consideration  or  in  press,  we'll  be  glad  to  send  page  proofs  with  a  note  of  the  embargo  date. 

If  your  press  office  wishes  to  issue  a  press  release,  they  should  contact  us  in  advance  for  final  embargo  information. 
Please  contact  our  press  office  at  press@cell.com.  Prior  to  publication,  please  discuss  any  possibility  of  the  paper 
being  referenced  elsewhere  in  the  literature  with  the  handling  editor. 

Again,  congratulations  on  a  very  nice  paper.  We  hope  you  found  the  review  process  to  be  constructive  and  are  pleased 
with  how  the  manuscript  was  handled  editorially.  We  look  forward  to  future  exciting  submissions  from  your  lab. 

Best  wishes, 

Sri 

Sri  Devi  Narasimhan,  Ph.D. 

Scientific  Editor,  Cell 
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23  Summary 

24  Despite  its  success  in  several  clinical  trials,  cancer  immunotherapy  remains  limited  by  the 

25  rarity  of  targetable  tumor-specific  antigens,  tumor-mediated  immune  suppression,  and  tox- 

26  icity  triggered  by  systemic  delivery  of  potent  immunomodulators.  Here  we  present  a  proof- 

27  of-concept  immunomodulatory  gene  circuit  platform  that  enables  tumor-specific  expres- 

28  sion  of  immunostimulators,  which  could  potentially  overcome  these  limitations.  Our  design 

29  comprised  de  novo  synthetic  cancer-specific  promoters  and,  to  enhance  specificity,  an 

30  RNA-based  AND  gate  that  generates  combinatorial  immunomodulatory  outputs  only  when 

31  both  promoters  are  mutually  active.  These  outputs  included  an  immunogenic  cell-surface 

32  protein,  a  cytokine,  a  chemokine,  and  a  checkpoint  inhibitor  antibody.  The  circuits  trig- 

33  gered  selective  T-cell-mediated  killing  of  cancer  cells  but  not  of  normal  cells  in  vitro.  In  in 

34  vivo  efficacy  assays,  lentiviral  circuit  delivery  mediated  significant  tumor  reduction  and 

35  prolonged  mouse  survival.  Our  design  could  be  adapted  to  drive  additional  immunomodu- 

36  lators,  sense  other  cancers,  and  potentially  to  treat  other  diseases  that  require  precise  im- 

37  munological  programming. 
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38  Introduction 


39  The  potential  of  cancer  immunotherapy  has  been  demonstrated  in  several  clinical  trials. 

40  Yet  this  approach  is  still  limited  by  several  major  challenges  (Maude  et  ah,  2014;  Postow  et  al., 

41  2015).  For  example,  CAR-T  cells  and  bispecific  antibodies  should  ideally  target  cell-surface  an- 

42  tigens  that  are  exclusively  present  on  tumor  cells,  since  targeting  normal  tissues  can  result  in  se- 

43  vere  side  effects  (Morgan  et  al.,  2010;  Sadelain  et  al.,  2013;  Satta  et  al.,  2013).  However,  finding 

44  highly  tumor-specific  surface  antigens  is  difficult,  which  limits  the  range  of  targetable  tumors 

45  (Hassan  and  Ho,  2008;  Kershaw  et  al.,  2014;  Klebanoff  et  al.,  2016;  Maude  et  al.,  2014;  Morgan 

46  et  al.,  2010).  In  addition,  even  when  ideal  targetable  tumor-antigens  are  available,  tumor- 

47  mediated  immunosuppression  can  prevent  successful  immunotherapy  by  disrupting  important 

48  immunological  functions  that  are  necessary  for  effective  anti-tumor  immune  responses 

49  (Chmielewski  et  al.,  2014;  Newtek  et  al.,  2016;  Rabinovich  et  al.,  2007;  Turley  et  al.,  2015). 

50  Thus,  the  lack  of  targetable  tumor- specific  antigens  and  tumor-mediated  immunosuppression 

5 1  constitute  major  challenges  for  cancer  immunotherapy. 

52  Immunostimulatory  factors  can  be  utilized  to  overcome  these  challenges,  but  off-target 

53  activity  often  results  in  severe  toxicity  (van  der  Burg  et  al.,  2016).  For  example.  Surface  T-cell 

54  Engagers  (STEs)  are  artificial  immunogenic  cell-surface  proteins  that  bind  the  non-variable  re- 

55  gions  of  the  T-cell  receptor  complex.  Consequently,  STE  expressing  cells  are  designated  to  T- 

56  cell  mediated  killing,  regardless  of  T-cell  receptor  antigen  specificity  (see  Figure  IB  for  our  STE 

57  design  based  on  membrane-anchored  anti-CD3  scEv).  Tumor-specific  STE  expression  can  sub- 

58  stitute  for  targetable  tumor  antigens,  but  must  be  constrained  to  cancer  cells  to  avoid  damage  to 

59  healthy  tissues  (Liao  et  al.,  2003;  Liao  et  al.,  2000;  Yang  et  al.,  2007).  Other  immunomodulators, 

60  such  as  chemokines,  cytokines  and  immune  checkpoint  inhibitors  can  assist  in  overcoming  tu- 

61  mor- mediated  immunosuppression,  but  have  caused  severe  side  effects  in  clinical  trials  (Lasek  et 

62  al.,  2014;  Leonard  et  al.,  1997;  Lin  et  al.,  2014;  Postow  et  al.,  2015).  Eurthermore,  individual 

63  immunomodulators  are  frequently  insufficient  on  their  own  to  yield  strong  anti-tumor  efficacy 

64  (Mahoney  et  al.,  2015;  Morrissey  et  al.,  2016).  Combinatorial  immunotherapy  can  provide  sig- 

65  nificantly  stronger  efficacy,  but  also  increase  the  risk  and  severity  of  adverse  effects  (Boutros  et 

66  al.,  2016;  Carson  et  al.,  2000;  Lee  et  al.,  2009;  Liao  et  al.,  2003;  Liu  et  al.,  2011;  Paul  et  al., 

67  2002).  Tumor-localized  release  of  immunomodulators  could  potentially  decrease  these  systemic 
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68  side  effects  and  improve  the  therapeutic  efficacy  (Neri  and  Sondel,  2016;  Shukla  and  Steinmetz, 

69  2016;  Tugues  et  ah,  2015),  but  is  challenging.  Thus,  although  various  strategies  are  being  ex- 

70  plored  to  overcome  this  problem,  further  development  is  required  (Bell  and  McFadden,  2014; 

71  Neri  and  Sondel,  2016;  Shukla  and  Steinmetz,  2016;  Tugues  et  ah,  2015;  Zaharoff  et  ah,  2010). 

72  Synthetic  biology  enables  the  creation  of  powerful  genetic  biological  tools  for  studying, 

73  diagnosing,  and  treating  disease  (Bacchus  et  ah,  2013;  Bojar  and  Fussenegger,  2016;  Fenno  et 

74  ah,  2014;  Kotula  et  ah,  2014;  Pardee  et  ah,  2016;  Schukur  and  Fussenegger,  2016;  Schukur  et  ah, 

75  2015;  Wei  et  ah,  2012;  Xie  et  ah,  2016).  For  example,  synthetic  tumor-targeting  gene  circuits 

76  have  been  previously  designed  to  target  tumor  cells  with  enhanced  specificity,  but  have  not  lev- 

77  eraged  the  anti-tumor  potential  of  the  immune  system,  which  potentially  limits  their  efficacy  (Liu 

78  et  ah,  2014;  Morel  et  ah,  2016;  Nissim  and  Bar-Ziv,  2010;  Xie  et  ah,  2011).  We  hypothesized 

79  that  gene  circuits  could  be  designed  to  achieve  tumor-specific  production  of  multiple  immuno- 

80  modulators  and  thus  overcome  major  obstacles  in  current  immunotherapies  (Figure  1  A). 

81  Here,  we  demonstrate  the  potential  of  gene  circuits  for  cancer  immunotherapy.  Tumor- 

82  targeting  gene  circuits  often  require  promoters  with  high  cell-state  specificity  to  function  as  can- 

83  cer  sensors,  which  can  be  hard  to  find  (Selvakumaran  et  ah,  2001;  Yun  et  ah,  2008).  Moreover, 

84  native  mammalian  promoters  can  span  several  kb  and  are  therefore  difficult  to  fit  into  the  limited 

85  encoding  capacity  of  many  viral  vectors.  Thus,  we  constructed  compact  synthetic  promoters  de 

86  novo  that  exhibited  enhanced  specificity  for  human  cancer  cells  compared  with  previously  de- 

87  scribed  endogenous  cancer-specific  promoters. 

88  Discriminating  cancer  cells  based  on  the  integration  of  multiple  cancer  markers  signifi- 

89  cantly  enhances  targeting  specificity  (Liu  et  ah,  2014;  Morel  et  ah,  2016;  Nissim  and  Bar-Ziv, 

90  2010;  Roybal  et  ah,  2016;  Xie  et  ah,  2011).  Thus,  to  enhance  the  circuit  selectivity  over  a  single 

91  synthetic  promoter,  we  designed  a  synthetic  RNA-based  circuit  that  integrates  the  activity  of  two 

92  synthetic  ovarian-cancer  specific  promoters  and  generates  output  proteins  only  when  both  pro- 

93  moters  are  decidedly  active,  analogous  to  a  Boolean  AND  gate.  Specifically,  input  promoter  1 

94  regulates  an  auto-inhibitory  mRNA  that  encodes  the  output  protein,  while  input  promoter  2  regu- 

95  lates  an  RNA  molecule  that  reliefs  the  auto-inhibition  and  thus  enables  output  production.  This 

96  RNA-based  design  has  the  potential  to  reduce  the  off-target  immunogenic  signature  of  the  circuit 

97  by  minimizing  the  expression  of  foreign  proteins  in  normal  cells.  Immunomodulatory  outputs 
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98  used  in  this  study  include  a  cell-surface  antigen  (STE)  together  with  secreted  a  chemokine  that 

99  promotes  T-cell  trafficking  to  tumor  sites  (CCL21),  a  cytokine  that  supports  T-cell  activation  and 

100  expansion  (IL12),  and  an  immune  checkpoint  inhibitor  (anti-PDl  antibody),  hereafter  referred  to 

101  as  the  SCIP  combination  (Lasek  et  ah,  2014;  Liao  et  ah,  2003;  Liao  et  al.,  2000;  Lin  et  al.,  2014; 

102  Postow  et  al.,  2015;  Yang  et  al.,  2007). 

103  We  characterized  our  circuits  with  in  vitro  specificity  assays  using  normal  and  cancer  cell 

104  lines  and  in  vivo  efficacy  assays  using  a  disseminated  ovarian  cancer  model,  which  recapitulates 

105  advanced  human  ovarian  cancer  (Lengyel  et  al.,  2014;  Sale  and  Orsulic,  2006).  In  addition,  we 

106  demonstrated  that  our  circuit  can  be  readily  modified  to  target  breast  cancer  cells  by  replacing 

107  the  synthetic  input  promoters.  Our  circuit  design  is  modular  and  enables  facile  construction,  test- 

108  ing,  and  optimization  of  arbitrary  promoter  inputs,  circuit  components,  and  outputs.  The  outputs 

109  can  be  any  genetically  encodable  element,  including  intracellular,  secreted,  and  cell-surface  pro- 

110  teins.  This  flexibility  will  enable  circuits  to  be  optimized  using  advanced  tumor  models  and  pa- 

111  tient-derived  tumors,  which  will  be  important  for  translation  into  clinical  use.  Thus,  our  pro- 

112  grammable  immunomodulatory  gene  circuits  have  broad  potential  utility  for  cancer  immunother- 

113  apy  and  studying  tumor  immunology. 

1 14  Results 

115  RNA-based  AND  Gate  Circuit 

116  We  first  designed  and  optimized  the  RNA-based  AND  gate.  We  engineered  two  genetic 

117  modules  such  that  each  was  regulated  by  a  separate  promoter  (Figure  2A-B).  Module  1  consists 

118  of  input  promoter  1  (“PI”),  which  regulates  the  expression  of  an  RNA  that  encodes  two  exons  of 

1 19  an  output  protein  mKate2  (mK-Exl,  mK-Ex2),  separated  by  a  synthetic  intronic  miRNA  (miRl), 

120  and  perfect  match  miRl  binding  sites  (BS(Pe))  downstream  to  the  second  exon.  Following  splic- 

121  ing,  the  exons  are  assembled  into  a  mature  output  mRNA  that  contains  BS(Pe)  in  the  3’  untrans- 

122  lated  region  (3’  UTR)  while  the  intron  is  processed  by  the  cell  to  produce  miRl.  Consequently, 

123  when  only  promoter  1  is  active,  miRl  binds  the  BS(Pe)  sites  in  the  output  mRNA  and  targets  it 

124  for  degradation.  Therefore,  Module  1  constitutes  an  auto-inhibitory  loop  for  the  output  protein 

125  production  (Figure  2A).  Module  2  consists  of  input  promoter  2  (“P2”),  which  regulates  a  miRNA 

126  sponge  that  consists  of  multiple  bulged  complementary  miRl  binding  sites  (BS(B))  and  thus  acts 
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127  as  a  miRl  shunt  that  relieves  the  auto-inhibitory  nature  of  Module  1  when  it  is  expressed  (Figure 

128  2A-B).  In  this  system,  the  output  protein  is  expressed  at  a  high  level  only  when  the  promoters 

129  regulating  both  modules  are  mutually  active.  However,  when  only  one  of  the  promoters  is  active 

130  or  when  none  of  the  promoters  are  active,  output  protein  expression  is  low. 

131  We  systematically  optimized  various  aspects  of  the  miRNA-based  mechanism  (Figure  3 

132  &  Figure  Sl-3;  for  further  details,  see  Supplementary  Information)  (Auyeung  et  al.,  2013;  Ebert 

133  and  Sharp,  2010;  Fellmann  et  al.,  2013;  Cure  et  al.,  1997;  Haraguchi  et  al.,  2009;  Kertesz  et  al., 

134  2007;  Rogakou  et  al.,  1998;  Suzuki  et  al.,  2015).  In  Module  2,  we  characterized  how  the  miRNA 

135  binding  site  sequence,  expression  level,  and  sponge  architecture  determines  the  efficiency  of  the 

136  miRl  sponge  (Figure  3  &  Figure  Sl-2).  We  defined  four  different  states  for  this  circuit:  state 

137  [0,0]  -  both  Module  1  and  Module  2  are  inactive,  [0,1]  -  Module  1  is  inactive  while  Module  2  is 

138  active,  [1,0]  -  Module  1  is  active  while  Module  2  is  inactive,  and  [1,1]  -  both  Module  1  and 

139  Module  2  are  active  (Figure  2B).  The  ON-OFF  ratio  for  the  AND  gate  circuit  was  defined  as  the 

140  mKate2  output  level  in  the  presence  versus  the  absence  of  the  sponge  (state  [1,1]  versus  state 

141  [1,0]  in  Figure  2B).  Ultimately,  we  chose  a  sponge  architecture  Sponge  7:  H2Alp-ECl- 

142  10Bzranbi-EC2-5B  (Figure  S2A)  as  the  final  ‘optimized  sponge’  design,  since  it  yielded  the 

143  highest  ON-OFF  ratio  for  the  AND  gate  circuit  (~6-fold.  Figure  S2A). 

144  In  Module  1,  we  optimized  the  synthetic  intronic  pri-miRNA  backbone  in  the  self- 

145  inhibitory  mKate2  transcript  for  enhanced  miRl  production  efficiency  (Greber  et  al.,  2008; 

146  Nissim  et  al.,  2014;  Xie  et  al.,  2011).  The  absolute  mKate2  signal  in  the  OFF  (no  sponge,  state 

147  [1,0])  and  ON  states  (optimized  sponge,  state  [1,1])  were  both  affected  by  the  pri-miRNA  de- 

148  signs  we  tested  (Mvl-Mv3,  Figure  S3).  Mvl  exhibited  an  OFF-ON  shift  from  -1400  to  -8400 

149  Mean  Fluorescence  Intensity  (MFI)  (5.8-fold  activation  ratio),  Mv2  from  -170  to  -500  MFI  (3- 

150  fold  activation),  and  Mv3  from  -380  to  -2600  MFI  (6.8-fold  activation).  We  selected  the  Mv3 

151  architecture  for  further  experimentation  since  it  had  low  mKate2  level  in  OFF  state,  substantial 

152  absolute  mKate2  levels  in  the  ON  state,  and  the  highest  fold-activation.  However,  we  note  that 

153  different  designs  could  be  more  suitable  for  other  targeting  scenarios.  For  example,  when  using 

154  effector  outputs  with  low  potency,  Mvl  could  be  chosen  to  achieve  high  absolute  levels  of  ex- 

155  pression  in  the  ON  state.  Conversely,  Mv2  could  be  suitable  for  expressing  extremely  potent 

156  outputs  in  which  even  low  background  output  activity  could  be  detrimental. 
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158  Synthetic  Cancer-Specific  Promoters 

159  To  generate  synthetic  ovarian-cancer-specific  promoters,  we  identified  transcription  fac- 

160  tors  (TFs)  that  are  overexpressed  in  ovarian  cancer  cells  by  comparing  gene  expression  in  prima- 

161  ry  human  ovarian  tumors  versus  healthy  controls,  and  their  respective  binding  motifs,  using  pub- 

162  lically  available  databases  (for  more  details,  see  Supplementary  Information).  In  each  synthetic 

163  promoter  S(TF)p,  multiple  binding  motifs  for  a  single  cancer-specific  TF  were  encoded  in  tan- 

164  dem  upstream  to  a  synthetic  minimal  promoter  derived  from  the  major  adenoviral  late  promoter 

165  (Figure  4A).  This  design  approach  resulted  in  compact  promoters  with  evenly  distributed  and 

166  dense  TF-BSs  corresponding  to  a  small  number  of  TF  families.  For  example,  analysis  of  the  1  kb 

167  upstream  of  the  transcriptional  start  site  (TSS)  of  the  native  tumor-specific  promoter  H2Alp 

168  (Rogakou  et  al.,  1998)  revealed  ~82  possible  binding  sites  for  40  different  TF  families,  with  an 

169  mean  of  ~8  TF-BSs  per  lOObp.  In  comparison,  the  synthetic  S(E2Fl)p  promoter,  comprised  of 

170  E2F1  binding  sites,  had  only  44  possible  binding  sites  for  7  different  TF  families,  with  a  mean  of 

171  ~16  TF-BSs  per  lOObp  (Table  SI;  for  more  details  on  promoters  analysis,  see  Supplementary 

172  Information). 

173  We  assayed  the  specificity  and  activity  of  natural  and  artificial  promoters  by  infecting 

174  human  OVCAR8  (OV8)  ovarian  cancer  cells  and  normal  primary  cells,  specifically  primary 

175  adult  human  dermal  fibroblasts  (aHDF)  and  primary  human  ovarian  epithelium  (HOV).  Several 

176  synthetic  promoters,  including  S(E2Fl)p  and  S(cMYC)p  (comprised  of  E2F1  and  cMYC  binding 

177  sites,  respectively),  were  highly  active  in  OVCAR8  ovarian  cancer  cells  but  not  in  normal  cells. 

178  These  synthetic  promoters  exhibited  greater  activation  ratios  compared  with  known  native  tu- 

179  mor-specific  promoters  such  as  SSXlp  (Gure  et  al.,  1997)  and  H2Alp  (Rogakou  et  al.,  1998) 

1 80  (Figure  4B  for  selected  promoters). 

181  For  example,  SSXlp  exhibited  9-fold  activation  in  OVCAR8  compared  to  aHDF  and  17- 

182  fold  activation  compared  to  HOV-epi  (Figure  4B).  H2Alp  had  2-fold  activation  in  OVCAR8 

183  compared  to  aHDF  and  5.6-fold  compared  to  HOV.  The  synthetic  promoter  S(E2Fl)p  had  224- 

184  fold  induction  in  OVCAR8  versus  aHDF  and  570-fold  activation  in  OVCAR8  compared  to 

185  HOV,  with  absolute  maximal  expression  levels  comparable  to  the  native  promoters.  S(cMYC)p 

186  exhibited  15-fold  induction  in  OVCAR8  versus  aHDF  and  36-fold  activation  in  OVCAR8  com- 
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187  pared  to  HOV.  Since  the  maximal  activity  of  S(cMYC)p  was  lower  than  S(E2Fl)p  in  OVCAR8 

188  cells,  we  selected  S(cMYC)p  to  drive  expression  of  Module  1  and  S(E2Fl)p  to  drive  the  expres- 

189  sion  of  Module  2.  This  was  done  so  that  the  sponge  could  be  expressed  at  a  greater  level  than  the 

190  miRl -encoding  construct,  which  is  important  for  achieving  a  high  activation  ratio  with  our  cir- 

191  cuit  (Figure  S2A). 

192 

193  Multi-Output  AND  Gate  for  Combinatorial  Immunomodulation 

194  To  facilitate  the  expression  of  combinatorial  immunomodulators  from  our  circuit,  we  re- 

195  placed  mKate2  in  Module  1  with  a  synthetic  transcription  factor  that  activates  the  expression  of 

196  downstream  genes  regulated  by  cognate  synthetic  promoters  (Figure  5A).  The  artificial  transcrip- 

197  tion  factor  implemented  was  the  yeast  GAL4  DNA-binding  domain  fused  to  the  viral  VP  16  tran- 

198  scriptional  activation  domain  (GAD),  although  other  transcription  factors  could  be  used  instead. 

199  Specifically,  we  used  the  synthetic  promoter  S(cMYC)p  to  express  Module  1  (S(cMYC)p- 

200  GADExi-[niiRl]cons-GADEX2-llPe).  Module  1  encoded  a  self-inhibitory  RNA  transcript  with 

201  miRl  as  an  intron  within  the  GAD  coding  sequence  followed  by  11  tandem  perfect  match  miRl 

202  binding  sites  (llPe)  in  the  3’  UTR.  In  the  Module  2  vector,  the  synthetic  promoter  S(E2Fl)p 

203  regulated  the  expression  of  the  optimized  sponge  architecture  (Sponge  7:  EC1-10Bzranbi-EC2- 

204  5B).  Finally,  we  encoded  an  mKate2  output  under  the  regulation  of  an  artificial  promoter  target- 

205  ed  by  GAD  that  contains  GAL4-BSs  upstream  to  a  synthetic  minimal  late  adenoviral  promoter 

206  (GALp),  namely  Module  3  (Figure  5A).  We  hypothesized  that  increasing  the  number  of  GAL4- 

207  BSs  in  the  GALp  would  provide  a  convenient  strategy  for  increasing  output  expression  levels. 

208  We  therefore  built  3  different  versions  of  GALp,  either  with  5,  8  or  14  tandem  GAL4-BSs  (G5p, 

209  G8p,  and  G14p,  respectively).  We  further  surmised  that  output  levels  in  the  OFF  state  could  be 

210  lowered  by  adding  miRl-BS(Pe)  to  the  mKate2  output,  thus  allowing  miRl  to  reduce  both  GAD 

21 1  and  mKate2  levels  directly  when  the  sponge  is  not  active  (Figure  5A).  We  therefore  built  3  addi- 

212  tional  versions  of  the  output  vector  in  which  a  single  perfect  match  miRl-BS  (IPe)  was  encoded 

213  in  the  3’  UTR  of  the  output  mKate2  transcript  (G5-Pe,  G8-Pe  and  G14-Pe,  respectively).  This 

214  design  architecture  enables  convenient  tuning  of  each  output’s  expression  levels  based  on  the 

215  promoter  and  whether  IPe  was  incorporated  into  the  output  transcript. 
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216  We  first  characterized  the  multi-output  AND-gate  architectures  by  modeling  scenarios 

217  when  neither  (state  [0,0]),  only  one  (states  [0,1]  or  [1,0]),  or  both  cancer-specific  promoters 

218  (state  [1,1])  are  active  in  a  given  cell  type  (Figure  5B).  To  do  so,  we  encoded  each  module  on 

219  separate  lentiviruses  and  infected  OVCAR8  cells  with  the  corresponding  Module  3  architecture 

220  and  either:  (A)  a  negative  control  Module  1  that  expresses  the  non-specific  transcription  factor 

221  rtTA3  {Module  Icon)  and  a  negative  control  Module  2  in  which  S(E2Fl)p  expresses  ECFP  with- 

222  out  miRl-BS(B)  {Module  Icon)  to  model  the  scenario  when  neither  cancer-specific  promoters 

223  are  active  [0,0];  (B)  Module  Icon  and  Module  2  to  model  the  scenario  when  only  S(E2Fl)p  is 

224  active  [0,1];  (C)  Module  1  and  Module  2con  to  model  the  scenario  when  only  S(cMYC)p  is  ac- 

225  tive  [1,0];  (D)  Module  1  and  Module  2  to  model  the  scenario  when  both  promoters  are  active 

226  [1,1], 

227  The  circuits  expressed  minimal  mKate2  fluorescence  in  states  [0,0]  and  [0,1].  In  state 

228  [1,0],  there  was  a  low  background  mKate2  level  that  increased  from  5  to  8  GAL4-BSs  in  the  syn- 

229  thetic  GALp  and  was  attenuated  with  the  presence  of  miRl-BS(Pe)  in  Module  3  (Eigure  5B).  On- 

230  ly  when  both  input  constructs  were  present  (state  [1,1])  did  we  observe  a  significantly  increased 

231  mKate2  output.  Furthermore,  the  mKate2  output  level  could  be  enhanced  by  increasing  the  num- 

232  her  of  GAL4-BSs  in  the  Module  3  synthetic  promoter  from  5  to  8  GAL4-BSs.  Although  adding 

233  perfect  match  miRl  binding  sites  in  the  3’  UTR  of  the  mKate2  transcript  decreased  the  basal  and 

234  maximal  output  levels,  it  also  improved  fold-activation  (G5-Pe,  G8-Pe,  and  G14-Pe  versus  G5, 

235  G8,  and  G14,  respectively).  The  G8-Pe  circuit  exhibited  negligible  mKate2  levels  in  states  [0,0] 

236  and  [0,1],  and  high  8-fold  mKate2  activation  in  state  [1,1]  over  the  minimal  mKate2  expression 

237  in  state  [1,0],  and  was  thus  chosen  for  further  use. 

238  We  next  compared  mKate2  expression  from  the  G8-Pe  circuit  delivered  via  lentivirus  in 

239  OVCAR8  cells  versus  normal  cells,  including  primary  aHDF,  HOV,  and  human  T  cells.  The  cir- 

240  cuit  exhibited  higher  mKate2  expression  in  OVCAR8  versus  aHDF  and  HOV  in  the  [1,1]  circuit 

241  configuration  (4-fold  and  71-fold,  respectively.  Figure  5C).  We  next  determined  the  activity  of 

242  the  S(cMYC)p  and  S(E2Fl)p  promoters  individually,  as  well  as  the  activity  of  the  G8-Pe  circuit 

243  with  S(cMYC)p  and  S(E2Fl)p  promoters  as  inputs,  in  primary  human  T  cells.  We  observed  that 

244  S(E2Fl)p  was  significantly  activated  in  primary  human  T  cells.  However,  since  S(cMYC)p  was 

245  not  active  in  these  cells,  the  integrated  output  generated  by  the  G8-Pe  circuit  was  30-fold  lower 
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246  compared  to  the  activity  of  S(E2Fl)p  on  its  own  (Figure  5D  &  S4A).  Furthermore,  we  tested  the 

247  specificity  of  this  design  using  additional  cell  lines,  including  non-tumorigenic  human  ovarian 

248  fibroblasts  (HOF),  non-tumorigenic  human  ovarian  microvasculature  epithelial  cells  (HOMFC), 

249  non-tumorigenic  colonic  epithelial  cells  (CCD-841 -CoN),  and  non-tumorigenic  mammary  epi- 

250  thelial  cells  (MCF-lOA  and  MCF-12A),  as  well  as  inducible  pluripotent  stem  cells  (iPSCs).  The 

251  circuit  achieved  high  specificity  for  OVCAR8  cells,  with  10-fold  to  800-fold  higher  circuit  out- 

252  put  in  OVCAR8  cells  versus  these  additional  cell  lines  (Figure  S4B).  These  data  demonstrate  the 

253  stringency  of  our  circuit  and  the  potential  of  our  circuit  design  to  enhance  therapeutic  specificity 

254  by  reducing  the  probability  of  false-positive  targeting  compared  to  single  promoter. 

255  Finally,  to  demonstrate  the  modularity  of  our  design,  we  examined  whether  the  G8-Pe 

256  circuit  could  be  adapted  to  target  breast  cancer  cells  by  simple  modification  of  the  input  promot- 

257  ers.  For  this  purpose,  we  created  two  additional  synthetic  promoters,  S(USFl)p  and  S(MAFK)p, 

258  that  were  designed  to  specifically  target  the  MDA-MB-453  breast  cancer  cell  line  but  not  the 

259  MCF-lOA  non-tumorigenic  mammary  epithelial  cell  line.  We  compared  the  mKate2  outputs 

260  generated  by  each  promoter  individually  versus  the  G8-Pe  circuit  in  both  cell  lines.  This  new  de- 

261  sign  generated  a  high  output  only  in  MDA-MB-453  cells,  demonstrating  the  versatility  of  our 

262  circuit  architecture  and  its  ability  to  be  adapted  to  a  variety  of  different  therapeutic  scenarios  and 

263  applications. 

264 

265  STEs  trigger  T-cell-mediated  cancer  cell  killing  and  IFN-y  secretion 

266  Next,  we  sought  to  encode  immunotherapies  as  the  outputs  of  our  synthetic  cancer- 

267  specific  circuits.  An  STE  was  constructed  by  fusing  an  scFv  targeting  human  CD3s  with  an  inert 

268  membrane  anchor  derived  from  the  Duffy  antigen/receptor  for  chemokines  (DARC)  of  erythro- 

269  cytes  (von  Nickisch-Rosenegk  et  ak,  2012)  (Figure  IB).  We  replaced  mKate2  with  STE  as  the 

270  output  of  the  G8-Pe  circuit  and  validated  the  therapeutic  specificity  of  this  circuit  in  vitro.  We 

271  observed  20  to  118-fold  higher  STE  expression  on  circuit-transduced  OVCAR8  cells  than  on 

272  aHDF  or  HOV  cells  (Figure  6A).  We  also  observed  robust  T-cell-mediated  killing  of  circuit- 

273  transduced  OVCAR8  cells,  which  was  significantly  higher  than  in  aHDF  or  HOV  cells  (Figure 

274  6B  &  Figure  S4C-E).  There  was  minimal  T-cell-mediated  killing  of  non-STE-expressing  cells 

275  (Figure  S4C-E),  which  could  be  caused  by  our  experimental  conditions,  including  MHC  mis- 
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276  match  between  T  cells  and  target  cells  as  well  as  T-cell  mediated  killing  with  non-STE-based 

277  mechanisms  (Kagi  et  ah,  1994).  In  addition  to  direct  tumor  killing,  robust  IFN-y  secretion  by  T 

278  cells  is  critical  for  strong  therapeutic  efficacy.  We  observed  strong  and  specific  IFN-y  secretion 

279  by  T  cells  only  when  they  were  co-cultured  with  circuit-transduced  OVCAR8  cells  (Figure  6C  & 

280  Figure  S4F),  which  was  20-28  fold  higher  than  when  T  cells  were  co-cultured  with  circuit- 

281  transduced  aHDF  or  HOV  cells.  We  also  observed  minimal  OVCAR8  killing  and  IFN-y  secre- 

282  tion  by  T  cells  when  Module  2  was  removed  from  the  G8-Pe  circuit  (state  [1,0])  to  model  the  po- 

283  tential  scenario  in  which  one  of  the  promoters  exhibits  some  off-target  activity  (Figure  6A-C  & 

284  Figure  S4C-F). 

285  We  also  performed  a  control  experiment  in  which  STE  expression  was  regulated  by  the 

286  human  Ubiquitin  C  (hUbC)  promoter,  which  is  constitutively  active  in  OVCAR8,  aHDF,  and 

287  HOV  cells.  We  examined  the  sensitivity  of  these  cells  to  STE-mediated  killing  and  their  ability 

288  to  trigger  IFN-y  production  by  T  cells  (Figure  S4G-I).  We  found  that  while  STE  expression  lev- 

289  els  in  OVCAR8  cells  were  ~3-fold  higher  than  in  aHDF  and  HOV  cells,  STE-mediated  killing 

290  was  comparable  in  aHDF  and  OVCAR8,  and  ~2-fold  higher  in  HOV.  In  addition,  all  three  STE- 

291  expressing  cell  lines  triggered  comparably  strong  IFN-y  production  by  T  cells.  These  results 

292  highlight  the  stringency  of  our  AND  gate,  which  mediates  significantly  higher  OVCAR8  speci- 

293  ficity  over  these  other  cell  lines. 

294  We  next  examined  whether  STE  expression  could  trigger  T-cell-mediated  anti-tumor  ac- 

295  tivity  in  vivo.  OVCAR8  cells  were  transduced  with  lentiviruses  encoding  a  doxcycline-inducible 

296  STE  expression  system  and  then  injected  i.p.  into  immunodeficient  NSG  mice  on  day  0.  Since 

297  NSG  mice  lack  T  cells,  human  T  cells  were  also  periodically  injected  i.p.  Specifically,  STE 

298  expression  was  induced  at  day  14  followed  by  human  T  cell  injections  at  day  16  and  31.  In  the 

299  absence  of  T  cells,  tumor  growth  was  similar  between  the  STE-induced  and  non-induced  groups. 

300  STE  expression  led  to  a  6-fold  decrease  in  tumor  burden  in  mice  treated  with  T  cells  compared  to 

301  non-treated  mice  (Figure  S4J-L),  demonstrating  that  the  anti-tumor  activity  of  STEs  was  depend- 

302  ent  on  T  cells.  However,  consistent  with  previous  studies  (Liao  et  al.,  2003;  Liu  et  al.,  2011;  Paul 

303  et  al.,  2002),  STE  expression  on  its  own  was  unable  to  completely  abrogate  tumor  growth,  even 

304  when  expressed  in  high  levels  in  nearly  all  the  tumor  cells  (Figure  S4L). 

305 
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306  Tumor-specific  expression  of  multiple  immunotherapeutic  outputs  triggers  effective  T- 

307  cell-dependent  tumor  killing 

308  Since  combination  therapy  has  shown  enhanced  therapeutic  efficacy  over  monotherapies 

309  (Lasek  et  ah,  2014;  Mahoney  et  al.,  2015;  Morrissey  et  al.,  2016;  Tugues  et  al.,  2015)  and  STE 

310  expression  alone  was  insufficient  to  achieve  strong  anti-tumor  effects  in  vivo,  we  engineered  our 

311  optimized  cancer- specific  gene  circuit  to  generate  multiple  immunomodulators  as  outputs.  To 

312  demonstrate  the  capability  of  the  circuit  to  generate  various  types  of  commonly  used  immuno- 

313  modulators,  we  implemented  the  SCIP  combination,  which  includes  a  cell-surface  antigen  (STE) 

314  together  with  secreted  chemokine  (CCL21),  cytokine  (IL12),  and  checkpoint  inhibitor  (anti-PDl 

315  antibody). 

316  First,  we  transduced  OVCAR8  cells  with  lentiviruses  encoding  the  G8-Pe  AND  gate  and 

317  the  SCIP  combination.  We  validated  the  selective  co-production  of  STE,  CCL21,  IL12,  and  an 

318  anti-PDl  antibody  from  OVCAR8  cells,  but  not  from  aHDF  or  HOV  cells,  in  vitro  (Figure  6D- 

319  G).  Next,  we  validated  the  functionality  of  our  SCIP  circuit  with  an  in  vivo  model  of  human 

320  ovarian  cancer.  We  engineered  the  G8-Pe  circuit  to  produce  a  control  output,  STE  only,  or  the 

321  SCIP  combination  in  OVCAR8  cells.  We  then  established  i.p.  disseminated  ovarian  cancers  by 

322  injecting  these  OVCAR8  cells  into  the  peritoneal  cavity  of  NSG  mice.  This  experimental  setup 

323  was  used  to  model  the  functionality  of  the  circuit  throughout  the  course  of  the  treatment  and  its 

324  effects  on  circuit  efficacy.  For  example,  the  activity  of  our  input  promoters  could  be  altered  in 

325  response  to  environmental  signals  or  due  to  tumor  heterogeneity,  resulting  in  circuit  inactivation 

326  and  reduced  efficacy. 

327  In  the  absence  of  T  cells,  all  tumors  grew  robustly  (Figure  S5A).  When  administered  with 

328  T  cells,  tumor  sizes  in  the  control  output  (Control)  and  STE  only  groups  were  significantly  larger 

329  than  the  STE-i-CCL21+IL12+anti-PDl  Ab  (SCIP)  combination,  in  which  tumors  were  undetecta- 

330  ble  in  most  mice  (Figure  7A  &  S5A-B).  Significantly  prolonged  mice  survival  was  also  observed 

331  in  the  SCIP  group  (Figure  7B).  Thus,  combinatorial  immunotherapies  expressed  from  within  tu- 

332  mor  cells  can  mediate  strong  anti-tumor  effects. 

333  Next,  we  sought  to  confirm  that  the  SCIP-expressing  circuit  could  mediate  robust  effica- 

334  cy  even  when  it  was  only  introduced  into  a  small  fraction  of  tumor  cells.  This  experiment  was 

335  used  to  model  a  common  problem  with  tumor-targeting  therapies  in  which  the  therapeutic  circuit 
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336  is  activated  only  in  a  fraction  of  the  tumor.  For  example,  this  phenomenon  can  be  mediated  by 

337  modulation  of  gene  expression,  cellular  heterogeneity  within  tumors,  or  partial  circuit  delivery. 

338  We  mixed  wild-type  OVCAR8  cells  with  OVCAR8  cells  eontaining  the  SCIP-expressing  circuit 

339  in  defined  ratios.  Even  when  only  15%  or  30%  of  the  overall  tumor  was  composed  of  SCIP- 

340  expressing  OVCAR8  eells,  signifieant  tumor  reduction  and  survival  was  still  achieved,  only  in 

341  the  presence  of  T  cells  (Figure  7C-D  &  Figure  S6A-B).  This  data  suggests  that  immunomodula- 

342  tory  gene  circuits  may  be  able  to  mediate  strong  anti-tumor  effects  despite  incomplete  delivery  in 

343  future  applications. 

344  Finally,  we  direetly  tested  the  efficaey  of  our  SCIP  circuit  within  in  vivo  delivery  settings. 

345  To  minimize  the  total  number  of  viruses  required  to  encode  the  entire  SCIP  circuit  and  thus  sim- 

346  plify  the  delivery  process,  we  eombined  the  CCL21  and  IL12  output  modules  into  a  single  lenti- 

347  viral  construe!  by  encoding  both  outputs,  separated  by  a  P2A  peptide  (Szymczak  et  al.,  2004), 

348  under  the  regulation  of  a  single  G8  output  promoter.  We  similarly  combined  the  anti-PDl  Ab 

349  heavy  chain  and  light  chain  modules  into  a  single  lentivirus.  Mice  with  intraperitoneally  dissem- 

350  inated  tumors  were  injected  i.p.  with  either  a  eontrol  lentivirus  mix  that  eontained  Module  1, 

351  Module  2  and  a  negative  control  output  or  a  therapeutic  virus  mix  that  contained  Module  1, 

352  Module  2  and  the  SCIP  outputs.  Human  T  cells  were  also  periodically  injected  i.p.  Under  these 

353  settings,  we  observed  significantly  reduced  tumor  burden  and  enhanced  survival  in  mice  treated 

354  with  the  SCIP  circuit  versus  the  control  (Figure  7E-F  &  Eigure  S7A-B).  This  data  indieates  that 

355  in  vivo  delivery  of  our  immunomodulatory  gene  cireuits  has  the  potential  to  trigger  anti-tumor 

356  effects,  even  when  it  is  encoded  on  multiple  vectors. 

357  Discussion 

358  In  this  study,  we  present  a  proof-of-concept  demonstration  that  versatile  genetic  circuits 

359  can  be  engineered  for  highly  speeifie  production  of  combinatorial  immunomodulators  from  with- 

360  in  cancer  eells.  This  modular  design  strategy  ean  be  further  adapted  to  target  additional  tumor 

361  types  and  optimized  for  clinical  translation.  For  example,  we  have  shown  that  the  circuit  can  be 

362  readily  modified  to  target  breast  cancer  cells  by  replacing  the  ovarian-cancer  specific  input  pro- 

363  moters  with  other  synthetic  promoters  that  are  breast-eancer  speeifie.  Furthermore,  we  showed 
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364  that  these  circuits  could  be  delivered  in  vivo  to  trigger  robust  anti-tumor  responses  and  increased 

365  survival  of  treated  mice. 

366  However,  the  in  vivo  efficacy  model  we  utilized  has  several  drawbacks  limiting  our  abil- 

367  ity  to  address  questions  such  as  in  vivo  targeting  specificity  and  safety  concerns.  In  our  tumor 

368  model,  tumor  cells  are  of  human  origin  and  all  other  normal  tissues  are  of  mouse  origin,  thus 

369  making  it  challenging  to  address  targeting  specificity.  Furthermore,  NSG  mice  injected  with  hu- 

370  man  T  cells  eventually  develop  life  threatening  graft-versus-host  disease  (GVHD)  (Covassin  et 

371  al.,  2011;  Schroeder  and  DiPersio,  2011),  making  it  difficult  to  study  safety  concerns  related  to 

372  therapy  in  this  model.  To  further  advance  this  approach  towards  therapeutic  applications,  these 

373  issues  must  be  addressed.  First,  design  and  testing  of  circuits  will  likely  need  to  be  performed  in 

374  additional  tumor  samples  and  tumor  models.  For  example,  specificity  and  efficacy  should  be  fur- 

375  ther  optimized  and  confirmed  using  patient-derived  cancer  cells  and  a  panel  of  normal  tissues. 

376  The  potential  off-target  activity  of  circuit  components,  such  as  miRl  and  the  miRNA  sponge, 

377  should  also  be  examined,  for  example  via  transcriptomics.  Moreover,  while  our  in  vitro  assays 

378  allow  for  testing  of  short-term  specificity,  future  efforts  should  measure  long-term  specificity 

379  and  side  effects  within  in  vivo  settings.  Second,  additional  in  vivo  cancer  models  that  account  for 

380  tumor  heterogeneity  should  be  used  in  efficacy  assays,  such  as  patient-derived  xenografts  (Morel 

381  et  al.,  2016).  Third,  we  evaluated  the  anti-tumor  efficacy  of  our  combinatorial  immunomodula- 

382  tors  in  NSG  mice  with  supplemented  human  T  cells,  which  does  not  capture  the  full  immune  re- 

383  sponse.  Future  work  will  be  needed  to  optimize  and  characterize  specific  immunomodulatory 

384  combinations  with  high  efficacy  and  safety  in  additional  cancer  models,  including  humanized 

385  mice  or  syngeneic  mouse  models  with  intact  immune  systems.  Such  systems  may  be  able  to 

386  model  the  efficacy  and  safety  effects  of  our  therapeutic  approach  on  other  arms  of  the  immune 

387  system  and  on  important  phenomena  such  as  epitope  spreading  (el-Shami  et  al.,  1999).  Mice  in 

388  such  systems  do  not  develop  GVHD  and  may  enable  us  to  address  safety  concerns  related  to  the 

389  therapy.  Fourth,  the  exact  combination  of  immunomodulatory  outputs  will  likely  need  to  be  op- 

390  timized  for  different  tumor  types.  Finally,  improved  delivery  strategies  will  be  needed  to  convey 

391  our  immunomodulatory  gene  circuits  into  tumors.  These  genetically  engineered  circuits  could  be 

392  delivered  into  tumors  using  non-viral  carriers,  such  as  HPAE  (Zhou  et  al.,  2016),  as  well  as  on- 

393  colytic  viruses,  such  as  HSV-1  (Andtbacka  et  al.,  2015;  Zhang  et  al.,  2013).  Given  that  only  15% 

394  of  tumor  cells  had  to  express  the  circuit  to  achieve  robust  anti-tumor  responses  in  our  model 
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395  (Figure  7C-D  &  Figure  S6)  and  that  our  circuit  could  mediate  anti-tumor  effects  when  delivered 

396  in  vivo  by  multiple  lentiviruses  (Figure  7E-F  &  Figure  S7),  we  hypothesize  that  delivery  effi- 

397  ciencies  will  not  need  to  be  perfect  to  achieve  a  therapeutic  effect.  However,  we  note  that  deliv- 

398  ery  requirements  will  likely  vary  in  other  tumor  models. 

399  Furthermore,  it  will  be  useful  to  test  whether  our  approach  can  be  used  in  combination 

400  with  other  therapeutic  strategies  to  enhance  efficacy.  For  example,  our  approach  may  be  useful  in 

401  improving  therapy  with  Tumor-Infiltrating  Lymphocytes  (TIL)  and  CAR-T  cells  by  tumor- 

402  specific  secretion  of  immunomodulators  that  attract  immune  cells  to  the  tumor  site  and  activate 

403  them  while  disrupting  the  immunosuppressive  microenvironment.  For  example,  TILs  and  CAR- 

404  T  cells  engineered  to  secrete  IL12  exhibited  stronger  efficacy  than  non-IL12  secreting  counter- 

405  parts  (Koneru  et  al.,  2015;  Pegram  et  al.,  2012;  Zhang  et  al.,  2015).  However,  IL12-secreting 

406  TILs  mediated  significant  systemic  toxicity  in  a  clinical  trial  (Zhang  et  al.,  2015),  probably  since 

407  T  cells  can  patrol  throughout  the  body  and  thus  produce  IL12  systemically.  Thus,  it  may  be  safer 

408  to  use  non-IL12-secreting  TILs  or  CAR-T  cells  together  with  circuit-triggered  IL12  secretion 

409  from  within  tumors. 

410  Additionally,  recent  studies  have  shown  that  many  tumors  express  immunogenic  neoanti- 

411  gens  that  can  be  targeted  by  the  host  immune  system  (Schumacher  and  Schreiber,  2015).  We 

412  speculate  that  using  our  strategy  to  trigger  tumor  killing  in  combination  with  neoantigen-based 

413  immunotherapies  could  be  a  fruitful  approach  to  investigate  (Carreno  et  al.,  2015;  Cohen  et  al., 

414  2015;  Gros  et  al.,  2016;  Gubin  et  al.,  2015;  Schumacher  and  Schreiber,  2015),  for  example,  by 

415  triggering  robust  T-cell-mediated  killing  of  cancer  cells  and  the  release  of  such  neoantigens.  Fur- 

416  thermore,  secreted  immunomodulators  could  disrupt  immunosuppressive  tumor  microenviron- 

417  ments  to  enhance  neoantigen  presentation  by  antigen-presenting  cells  and  further  improve  anti- 

418  tumor  immune  responses. 

419  Since  our  approach  enables  tumor-specific  secretion  of  defined  combinations  of  im- 

420  munomodulators,  it  could  be  a  useful  framework  to  study  the  immunological  mechanisms  that 

421  underlie  tumor  biology.  Moreover,  synthetic  gene  circuits  have  been  shown  to  have  the  potential 

422  to  treat  immunological  disorders  by  sensing  extracellular  disease  markers  (Schukur  et  al.,  2015). 

423  Similarly,  our  approach  for  sensing  internal  cell  states  with  multiple  artificial  promoters  could  be 

424  adapted  to  treat  other  complex  diseases  that  require  highly  specific  and  multifactorial  program- 
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425  ming  of  immunological  functions  with  simple  modifications  to  the  synthetic  promoter  inputs  and 

426  immunomodulatory  outputs. 

427  In  summary,  modulating  the  immune  system  is  a  promising  approach  to  treating  complex 

428  diseases.  However,  multiple  disease-mediated  mechanisms  can  limit  the  effectiveness  of  mono- 

429  therapies  (Mahoney  et  al.,  2015;  Morrissey  et  al.,  2016).  Combinatorial  immunomodulators  can 

430  be  used  to  overcome  this  difficulty,  but  given  the  broad  effects  of  the  immune  system  on  human 

43 1  physiology,  it  will  be  important  to  constrain  the  activity  of  immunomodulators  to  specific  disease 

432  compartments  to  avoid  undesirable  side  effects  and  maximize  efficacy  (Boutros  et  al.,  2016; 

433  Carson  et  al.,  2000).  Thus,  we  envision  that  synthetic  gene  circuit  immunotherapies  that  are  high- 

434  ly  specific,  tunable,  and  amenable  to  combinatorial  effector  expression  will  provide  powerful 

435  strategies  for  engineering  the  immune  system  to  understand  and  treat  disease. 

436 
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Methods 


Plasmid  construction 

The  various  circuit  modules,  synthetic  promoters,  and  immunomodulatory  constructs 
were  built  using  conventional  restriction  enzyme  cloning  and  Gibson  assembly.  Plasmids  com¬ 
positions  used  in  each  experiment  are  detailed  in  supplemental  file  “Plasmids  by  figures.pdf’. 
All  DNA  sequences  used  in  this  study  are  detailed  in  GenBank  format  in  the  single  text  file  “All 
sequences. docx”.  Additionally,  individual  GenBank  and  SnapGene  sequence  files  for  the  plas¬ 
mids  in  this  study  can  be  downloaded  from  a  public  Dropbox  folder  using  the  link  below: 
https://www.dropbox.com/sh/wbq2v31x87vzppe/ AAA  vFM7g6blK35AvlzSreoGDa?dl=0 

Cell  culture  and  cell  lines 

HEK-293T,  adult  Human  Dermal  Fibroblast  (aHDF),  MDA-MB-453,  MCF-lOA,  MCF- 
12A,  and  CCD-841-CoN  cells  were  obtained  from  the  American  Type  Culture  Collection,  Rock¬ 
ville,  MD  (HFK293T,  Catalog  #CRL-3216;  aHDF,  Catalog  #PCS-201-012;  MDA-MB-453,  Cat¬ 
alog  #HTB-131;  MCF-lOA,  Catalog  #CRL-10317;  MCF-12A,  Catalog  #CRL-10782;  CCD-841- 
CoN,  Catalog  #CRL-1790).  Ovarian  carcinoma  cell  line  8  (OVCAR8)  engineered  to  stably  ex¬ 
press  firefly  luciferase  (OVCAR8-Luc)  was  a  gift  from  S.  N.  Bhatia  (Massachusetts  Institute  of 
Technology,  Cambridge,  MA).  Primary  human  ovarian  epithelial  cells  (HOV),  primary  human 
ovarian  fibroblasts  (HOF),  and  primary  human  ovarian  micro  vasculature  endothelial  cells 
(HOMEC)  were  obtained  from  ScienCell  Research  Laboratories,  Carlsbad,  CA  (HOV,  Catalog 
#7310;  HOF,  Catalog  #7330;  HOMEC,  Catalog  #7300).  Inducible  pluripotent  stem  cell  (iPSC) 
was  a  gift  from  R.  Weiss  (Massachusetts  Institute  of  Technology).  HEK-293T,  OVCAR8-Luc, 
aHDF,  and  MDA-MB-453  cells  were  cultured  in  DMEM  (Life  Technologies,  Carlsbad,  CA) 
supplemented  with  10%  fetal  bovine  serum  (FBS;  VWR,  Radnor,  PA;  Catalog  #95042-108),  1% 
Non-Essential  Amino  Acids  (MEM/NEAA;  Hyclone;  Catalog  #16777-186),  and  1%  Pen/Strep 
(Life  Technologies  Catalog  #15140-122)  at  37  °C  with  5%  C02.  HOV  cells  were  cultured  in 
Ovarian  Epithelial  Cell  Medium  (ScienCell  Research  Laboratories  Catalog  #7311)  supplemented 
with  Pen/Strep  (ScienCell  Research  Laboratories  Catalog  #0503),  Ovarian  Epithelial  Cell 
Growth  Supplement  (ScienCell  Research  Laboratories  Catalog  #7352)  and  grown  on  a  poly-L- 
lysine  (Sigma- Aldrich,  St.  Louis,  MO;  Catalog  #P4707-50ML)  coated-plate.  HOF  cells  were  cul- 
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468  tured  in  Fibroblast  Medium  (ScienCell  Research  Laboratories  Catalog  #2301)  supplemented  with 

469  Pen/Strep  (ScienCell  Research  Laboratories  Catalog  #0503),  Fetal  Bovine  Serum  (ScienCell  Re- 

470  search  Laboratories  Catalog  #0010),  Fibroblast  Growth  Supplement  (ScienCell  Re-search  Labor- 

471  atories  Catalog  #2352)  and  grown  on  a  poly-L-lysine  (Sigma- Aldrich  Catalog  #P4707-50ML) 

472  coated-plate.  HOMEC  cells  were  cultured  in  Endothelial  Cell  Medium  (ScienCell  Research  La- 

473  boratories  Catalog  #1001)  supplemented  with  Pen/Strep  (ScienCell  Research  Laboratories  Cata- 

474  log  #0503),  Fetal  Bovine  Serum  (ScienCell  Research  Laboratories  Catalog  #0025),  Endothelial 

475  Cell  Growth  Supplement  (ScienCell  Research  Laboratories  Catalog  #1052)  and  grown  on  a  Fi- 

476  bronectin  (Sigma- Aldrich  Catalog  #F1141-2MG)  coated-plate.  MCF-lOA  and  MCF-12A  cells 

477  were  cultured  in  MEGM  BulletKit  (Lonza,  Walkersville,  MD;  Catalog  #CC-3151  &  CC-4136). 

478  All  cell  lines  were  banked  directly  after  being  purchased  from  vendors  and  used  at  low  passage 

479  numbers.  Cell  lines  were  tested  negative  for  mycoplasma  contamination  with  IMPACT  PCR  by 

480  the  Division  of  Comparative  Medicine  Diagnostic  Lab  at  MIT. 

481  Human  PBMCs  were  obtained  from  leukoreduction  collar  (Brigham  and  Women’s  hospi- 

482  tal  Crimson  Core  Laboratory,  Boston,  MA)  with  gradient  centrifugation.  To  culture  and  expand 

483  human  T  cells,  human  PBMCs  were  stimulated  with  40  ng/mL  of  anti-human  CD3  antibody 

484  (clone:  OKT3;  BioLegend,  San  Diego,  CA  ;Catalog  #317304)  for  3  days  in  RPMI  1640  (Invitro- 

485  gen,  Carlsbad,  CA;  Catalog  #11875119),  supplemented  with  10%  heat-inactivated  FBS,  10  mM 

486  HEPES  (Life  Technologies  Catalog  #15630080),  0.1  mM  non-essential  amino  acids,  1  mM  sodi- 

487  um  pyruvate  (Life  Technologies  Catalog  #11360-070),  1%  Pen/Strep,  50  pM  2-Mercaptoethanol 

488  (Sigma- Aldrich  Catalog  #M3148-25ML  ),  and  50  lU/mL  recombinant  human  IL-2  (NCI,  Freder- 

489  ick,  MD). 

490 

491  Transfection  for  AND  gate  optimization 

492  To  optimize  the  AND  gate  circuit  in  Figure  2-3  and  Figure  Sl-3,  we  tested  different 

493  combinations  of  two  plasmids:  Module  1  variants  expressed  an  auto-inhibitory  mKate2  RNA 

494  transcript  regulated  by  synthetic  miRNA  (miRl)  and  Module  2  variants  expressed  sponges  for 

495  miRl  within  an  ECFP  transcript,  where  ECFP  enables  us  to  assess  sponge  expression.  In  brief, 

496  12  pi  of  FuGENE  HD  transfection  reagent  (Promega,  Madison,  WI;  Catalog  #E231 1)  mixed  with 

497  100  pi  of  OptiMEM  medium  (Life  Technologies  Catalog  #31985)  was  added  to  the  mixture  of 

498  two  plasmids  (1  pg  each).  During  20  minutes  incubation  of  FuGENE  HD/DNA  complexes  at 
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499  room  temperature,  HEK-293T  suspension  cells  were  prepared  and  diluted  to  1.2  x  10^  cells/ml  in 

500  culture  medium.  0.5  ml  of  diluted  cells  (0.6  x  10^  cells)  were  added  to  each  FuGENE  HD/DNA 

501  complex  tube,  mixed  well,  and  after  5  min  incubation  at  room  temperature  transferred  to  a  desig- 

502  nated  well  in  6-well  plate  containing  2  ml  cell  culture  medium,  followed  by  incubation  at  37  °C 

503  with  5%  C02.  Transfected  cells  were  prepared  for  FACS  analysis  at  48-hour  post-transfection. 

504 

505  Lentivirus  production  and  transduction 

506  Lentiviruses  were  produced  in  HEK-293T  cells  using  co-transfection  in  6-well  plate  for- 

507  mat.  In  brief,  12  pi  of  FuGENE  HD  was  mixed  with  100  pi  of  OptiMEM  medium  and  was  added 

508  to  a  mixture  of  3  plasmids:  0.5  pg  of  pCMV-VSV-G  vector,  0.5  pg  of  lentiviral  packaging 

509  psPAX2  vector,  and  1  pg  of  lentiviral  expression  vector.  During  20  minutes  incubation  of 

510  FuGENE  HD/DNA  complexes  at  room  temperature,  HEK-293T  suspension  cells  were  prepared 

511  and  diluted  to  3.6  x  10^  cells/ml  in  cell  culture  medium.  0.5  ml  of  diluted  cells  (1.8  x  10^  cells) 

512  were  added  to  each  FuGENE  HD/DNA  complex  tube,  mixed  well,  and  incubated  for  5  min  at 

513  room  temperature  before  being  added  to  a  designated  well  in  6- well  plate  containing  1  ml  cell 

514  culture  medium,  followed  by  incubation  at  37°C  with  5%  C02.  Media  of  transfected  cells  were 

515  replaced  with  2.5  ml  fresh  culture  media  18  hours  post  transfection.  Supernatant  containing  new- 

516  ly  produced  viruses  was  collected  at  48-hour  post-transfection,  and  filtered  through  a  0.45  pm 

517  syringe  filter  (Pall  Corporation,  Ann  Arbor,  MI;  Catalog  #4614). 

518  For  testing  of  synthetic  promoters  (Figure  4A-B),  1:9  dilutions  of  filtered  viral  superna- 

519  tants  were  prepared  and  3  ml  of  each  individual  virus  was  used  to  infect  2.5  x  10^  OVCAR8-Luc, 

520  aHDF  or  HOV  cells  in  the  presence  of  8  pg/ml  polybrene  (Sigma)  overnight.  For  optimizing  and 

521  testing  the  GAD-expressing  gate  (Figure  5A-E),  we  transduced  three  viruses  encoding  Module  1, 

522  Module  2,  and  Module  3.  Each  virus  was  diluted  1:3  and  we  used  3  ml  of  the  pooled  viruses  to 

523  infect  2.5  x  10^  target  cells  in  the  presence  of  8  pg/ml  polybrene  (Sigma)  overnight.  Cell  culture 

524  medium  was  replaced  the  next  day  after  infection  and  cells  were  cultured  for  at  least  one  week 

525  prior  to  FACS  analysis  or  cytotoxicity  assay.  For  creating  cell  lines  bearing  circuits  expressing 

526  multiple  outputs  (Figure  6D-G,  Figure  7),  two  rounds  of  lentiviral  infection  were  performed.  The 

527  viruses  for  the  first  round  consisted  of  Module  1,  Module  2,  and  an  output  construct  expressing 

528  STE.  The  viruses  for  the  second  round  consisted  of  output  constructs  expressing  all  other  out- 

529  puts. 
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530 

531  Flow  cytometry 

532  For  characterizing  fluorescent  protein  expression,  cells  were  resuspended  with  DMEM 

533  and  analyzed  by  a  LSRII  Fortessa  cytometer  (BD  Biosciences,  San  Jose,  CA).  To  characterize 

534  STE  display  and  expression  level  on  cell  surfaces,  various  cell  lines  were  stained  with  phyco- 

535  erythrin  labeled  anti-HA  tag  Ah  (Miltenyi  Biotec;  catalog  #130-092-257).  Data  analysis  was  per- 

536  formed  by  ElowJo  software  (TreeStar  Inc,  Ashland,  OR). 

537 

538  T  cell-mediated  cell  lysis  in  vitro 

539  T  cell-mediated  cell  lysis  was  measured  by  LDH  release  assay  with  CytoTox  96®  Non- 

540  Radioactive  Cytotoxicity  Assay  Kit  (Promega,  Madison,  WI;  catalog  #G1780).  Human  T  cells 

541  were  incubated  with  circuit-expressing  target  cells  at  a  effector  to  target  ratio  (E:T  ratio)  of  20:1 

542  (6x10^:3x10"^)  in  duplicate  wells  in  96  well  flat-bottom  plates.  T  cells  and  target  cells  were  also 

543  plated  alone  to  determine  spontaneous  LDH  release  from  T  cells  (“effector  spontaneous”)  and 

544  from  target  cells  (“target  spontaneous”).  A  set  of  target  cells  were  also  plated  and  lysed  with  ly- 

545  sis  buffer  for  determining  maximal  LDH  release  (“target  maximum”).  8.5  hours  later,  cell-free 

546  culture  supernatant  was  harvested  and  LDH  release  was  determined  following  the  manufactur- 

547  er’s  protocol.  T  cell-mediated  lysis  was  calculated  as:  percentage  of  T  cell-mediated  lysis  =  [(ex- 

548  perimental  O.D.  value  -  effector  spontaneous  O.D.  value  -  target  spontaneous  O.D.  value)  /  (tar- 

549  get  maximum  O.D.  value  -  target  spontaneous  O.D.  value)]  x  100.  Media  O.D.  values  were  sub- 

550  traded  from  all  samples  before  calculation.  STE  mediated  lysis  was  calculated  as:  T  cell  mediat- 

55 1  ed  lysis  of  STE  expressing  cells  -  T  cell  mediated  lysis  of  non-STE  expressing  cells.  The  results 

552  of  each  experiment  represent  three  biological  replicates. 

553 

554  IFN-y  production  by  T  cells 

555  Human  T  cells  were  co-cultured  with  tumor  cells  at  an  E:T  ratio  =  20:1  (6x10^:3x10"^)  in 

556  duplicate  wells  in  96  well  flat-bottom  plates.  Cell-free  medium  was  collected  after  24h  and  IFN- 

557  y  concentration  determined  by  Human  IFN-y  DuoSet  ELISA  (R&D  systems,  Minneapolis,  MN; 

558  catalog  #DY285). 

559 
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560  Multiple-output  circuit  triggers  STE,  CCL21,  IL12,  and  anti-PDl  Ab  production  spe- 

561  cifically  by  OVCAR8  cells 

562  To  validate  STE,  CCL21,  IL12,  and  anti-PDl  Ab  production  by  OVCAR8  cells  but  not 

563  by  normal  cells  (Figure  6D-G),  2x10^  OVCAR8-Luc,  HOV,  or  aHDF  cells  transduced  with  the 

564  SCIP-expressing  circuit  were  plated  in  6-well  plates  and  cell-free  medium  was  collected  after 

565  72h.  Displayed  STE  on  cell  surfaces  was  measured  via  staining  with  phycoerythrin  labeled  anti- 

566  HA  tag  Ab  (Miltenyi  Biotec)  or  phycoerythrin  labeled  anti-His  tag  Ab  (Miltenyi  Biotec)  and 

567  quantified  by  flow  cytometry  analysis.  CCL21  concentration  was  determined  by  Human 

568  CCL21/6Ckine  DuoSet  ELISA  (R&D  systems;  catalog  #DY366).  IL-12  concentration  was  de- 

569  termined  by  Human  IL-12  p70  DuoSet  ELISA  (R&D  systems;  catalog  #DY1270),  and  anti-PDl 

570  Ab  production  was  validated  by  a  competitive  staining  approach.  The  competitive  staining  was 

571  done  by  incubating  cell-free  medium  (circuit-conditioned  or  control  medium)  with  activated  hu- 

572  man  T  cells  for  30  minutes,  then  followed  by  staining  with  anti-human  CD279  (PD-1)  Ab  [clone: 

573  EH12.2H7;  BioLegend;  catalog  #329920].  The  percentage  signal  reduction  was  calculated  by  the 

574  following  equation:  (1  -  PDl  staining  MFI  of  T  cells  incubated  with  circuit-conditioned  medium 

575  /  PDl  staining  MFI  of  T  cells  incubated  with  control  medium)  *  100%. 

576 

577  Mice 

578  Female  NSG  mice  (NOD.Cg-PrMc*"'*^  //2rg''"^”'^VSzJ;  Stock  No:  005557)  were  purchased  from 

579  the  Jackson  Laboratory  (Bar  Harbor,  ME).  Mice  were  used  in  experiments  at  the  age  of  8-13 

580  weeks  old.  All  experiments  were  conducted  at  the  Koch  Institute  Animal  Facility  under  a  proto- 

581  col  approved  by  MIT’ s  Committee  on  Animal  Care. 

582 

583  In  vivo  experiments  validating  circuit  therapeutic  efficacy 

584  For  in  vivo  experiments,  animals  were  randomly  assigned  to  each  experimental  group.  No 

585  specific  inclusion/exclusion  criteria  were  used  for  animal  studies.  Investigators  were  not  blinded 

586  to  animal  groups.  Four  to  five  mice  were  used  in  each  experimental  group.  To  validate  STE  ther- 

587  apeutic  efficacy  on  ovarian  cancer,  OVCAR8-Luc  cells  expressing  doxycycline  (Dox)-inducible 

588  STE  (5x10^  cells)  were  injected  intraperitoneally  (i.p.)  into  NSG  mice  on  day  0,  and  STE  expres- 

589  sion  were  induced  by  Dox  [0.2mg/mL  Dox  +  2%  Sucrose  (Sigma- Aldrich;  catalog  #S3929-1KG) 

590  in  drinking  water]  at  day  14.  Activated  human  T  cells  (10^  cells)  were  injected  at  day  16  and  31. 
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591  For  measuring  tumor  burden,  mice  were  injected  with  30mg  of  VivoGlo™  Luciferin,  In  Vivo 

592  Grade  (Promega;  catalog  #P1043)  10  minutes  prior  to  bioluminescence  imaging  (IVIS  spectrum 

593  optical  imaging  system,  Xenogen,  Alameda,  CA). 

594  To  measure  the  efficacy  of  different  therapeutic  output  combinations  on  ovarian  cancer, 

595  OVCAR8-Luc  cells  (5x10^  cells)  engineered  to  express  various  outputs  were  injected  i.p.  into 

596  NSG  mice  on  day  0.  Activated  human  T  cells  (lO’  cells)  were  injected  at  day  7,  14,  21  and  28. 

597  Tumor  burden  was  monitored  by  IVIS  imaging  using  the  protocol  described  above.  To  measure 

598  survival,  mice  were  euthanized  when  they  developed  signs  of  distress,  such  as  ruffled  fur,  poor 

599  body  posture,  distended  abdomen,  and  jaundice. 

600  To  identify  the  minimal  percentage  of  tumor  cells  needed  to  be  transduced  with  the  SCIP- 

601  expressing  circuit  to  mediate  robust  therapeutic  efficacy,  OVCAR8-Luc  cells  engineered  with 

602  the  SCIP-expressing  circuit  were  mixed  with  wild- type  OVCAR8-Luc  cells  at  various  ratios 

603  (15%  to  85%;  30%  to  70%).  5x10^  mixed  OVCAR8-Luc  cells  were  injected  i.p.  into  NSG  mice 

604  on  day  0.  Activated  human  T  cells  (lO’  cells)  were  injected  at  day  9,  16,  and  23.  Tumor  burden 

605  was  monitored  by  IVIS  imaging  using  the  protocol  described  above.  We  used  the  criteria  de- 

606  scribed  above  to  measure  survival. 

607  To  deliver  the  SCIP  circuit  with  lentivirus  in  vivo,  lentiviruses  for  each  circuit  component 

608  (Module  1,  Module  2,  STE,  CCL21-I-IL12,  and  anti-PDl  Ab,  see  Figure  7E  for  more  details) 

609  were  produced  as  described  above.  For  the  SCIP  therapy  group,  10  mL  of  viruses  for  each  of  the 

610  five  circuit  components  were  collected  and  concentrated  15-fold  using  Amicon  lOOK  MW  cutoff 

611  ultrafiltration  columns  (Millipore,  catalog  #UFC910024).  For  the  control  group,  10  mL  of  virus- 

612  es  for  Module  i,  10  mL  of  viruses  for  Module  2,  and  30  mL  of  viruses  encoding  the  control  out- 

613  put  (rtTA3)  was  collected  and  concentrated.  After  concentration,  viruses  were  pooled  and  a  p24 

614  ELISA  kit  (Advanced  Bioscience  Laboratories,  catalog  #5421)  was  used  to  determine  the  viral 

615  titer.  The  pooled  viruses  had  titers  of  -7x10^*'  viral  particles/mL.  Polybrene  was  added  to  the  vi- 

616  ruses  to  reach  a  final  concentration  of  16  pg/mL.  600  pL  of  viral  supernatant  was  injected  into 

617  mice  intraperitoneally. 

618 
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Figure  1.  Immunomodulatory  synthetic  gene  circuits  are  programmed  to  be  selectively  activated  in  cancer  cells. 

(A)  Panel  1:  Immunomodulatory  synthetic  gene  circuits  are  designed  to  integrate  the  activity  of  two  tumor- 
specific  synthetic  promoters  (PI  and  P2)  with  an  RNA-based  AND  gate  mechanism  and  generate  combinatorial 
immunomodulator  outputs  only  when  input  promoters  are  mutually  active.  When  activated,  the  AND  gate 
expresses  a  synthetic  transcription  factor  (GAD:  a  fusion  protein  consisting  of  the  GAL4  DNA  binding  domain 
and  VP16  transcription  activating  domain),  which  drives  the  co-expression  of  combinatorial  immunomodulators, 
including  Surface  T-cell  Engagers  (STEs,  which  are  anti-CD3s  scFvs  displayed  on  the  cell  surface)  and  secreted 
CCL21,  IL12,  and  an  anti-PDl  antibody.  Panel  2:  The  circuits  are  triggered  to  express  immunomodulators  in 
cancer  cells  (red),  but  not  normal  cells  (blue)  (Black-line  circuit  diagrams  indicate  circuits  are  active  and  Grey¬ 
line  circuit  diagrams  indicate  circuits  are  inactive).  Panel  3:  Cancer-specific  expression  of  combinatorial 
immunomodulators  triggers  effective  T-cell  mediated  killing  of  the  cancer  cells. 

(B)  Schematic  drawing  of  STE  displayed  on  cell  surfaces.  DARC  refers  to  the  silent  Duffy  antigen/receptor  for 
chemokines. 
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Figure  2.  The  basic  RNA-only  single-output  AND  gate  design. 

(A)  Module  1  of  the  AND  gate  is  designed  as  an  auto-inhibitory  loop  such  that  it  represses  its  own  output. 
Module  2  is  designed  to  inhibit  the  auto-inhibition  of  Module  1.  Module  1  and  Module  2  are  regulated  by 
cancer- specific  promoters  PI  and  P2,  respectively.  The  output  from  Module  1  is  expressed  at  a  high  level  only 
when  both  PI  and  P2  are  active,  which  enhances  the  tumor  specificity  of  the  circuit. 

(B)  All  4  possible  input  states  and  their  respective  output  states  are  shown  for  the  RNA-only  single-output  AND 
gate.  Input  states  are  defined  within  the  square  brackets  by  whether  Module  1  and  Module  2  are  active,  where  0 
means  inactive  and  1  means  active.  In  State  [1,0],  PI  is  active  and  the  mKate2  transcript  is  expressed  (mK-Exl 
and  mK-Ex2  denote  mKate2  exonl  and  mKate2  exon2,  respectively).  However,  this  transcript  encodes  an 
intronic  synthetic  miRNA  (miRI)  within  the  mKate2  gene  that  inhibits  the  mKate2  transcript  by  targeting 
perfect  match  miRI  binding  sites  (BS(Pe))  in  the  3’  untranslated  region.  Thus,  the  mKate2  levels  are  minimal. 
In  State  [0,1],  PI  is  inactive  so  the  output  protein  mKate2  is  not  expressed.  In  State  [0,0],  neither  PI  nor  P2  are 
active  and  thus  the  output  protein  mKate2  is  not  expressed.  In  State  [1,1],  the  P2  promoter  expresses  a  sponge 
for  miRI  that  is  based  on  bulged  miRI  binding  sites  (BS(B)).  This  enables  sequestration  of  miRI  away  from 
inhibiting  the  mKate2  transcript  expressed  by  PI,  thus  allowing  for  mKate2  expression. 


Figure  3 


Figure  3.  Flowchart  of  circuit  optimization  strategies  for  the  RNA-based  AND  gate. 

Specific  details  are  described  in  Supplementary  Information  and  Figure  SI -3.  Multiple  versions  of  the  miRl 
production  backbone  in  Module  1  as  well  as  the  miRl  binding-site  sequences  (BSs)  and  sponge  architecture  in 
Module  2  were  systematically  tested.  The  optimized  circuit  utilized  the  miR-E  backbone  (Mv3,  Table  S2)  in 
Module  1  to  generate  miRl  and  a  sponge  architecture  in  Module  2  that  encoded  10  optimized  bulged  miRl 
binding  sites  within  an  ECFP  gene  (EC1-10Bzranbi-EC2)  and  5  additional  bulged  miRl  binding  sites  (5B)  in 
the  downstream  3’UTR  (Figure  S2A).  The  mKate2  output  levels  for  the  optimized  circuit  with  a  control  sponge 
with  no  miRl  binding  sites  (‘Mv3  +  Cont.’,  representing  the  [1,0]  state  in  Figure  2C),  and  the  optimized  circuit 
with  the  optimized  sponge  (‘Mv3  +  Sponge’,  corresponding  to  Sponge  7  from  Figure  S2A,  representing  the  [1,1] 
state  in  Figure  2C)  are  shown.  Module  1  was  expressed  by  SSXlp  and  Module  2  was  expressed  by  H2Alp. 
Error  bars  represent  S.E.M.,  N  =  3  biological  replicates.  MFI  =  Mean  Fluorescence  Intensity.  BS  =  binding  sites. 
TuD  =  tough-decoy  architecture. 


See  also  Figure  SI -3. 


Figure  4 


A 


Synthetic  promoter  S(TF)p 


TF-BS  TF-BS  TF-BS  TF-BS 


lateADEp 


B 


Figure  4.  Synthetic  tumor-specific  promoters  exhibited  high  tumor  specificity. 

(A)  Design  of  synthetic  tumor-specific  promoters,  S(TF)p.  Multiple  binding  sites  for  a  transcription  factor  that 
is  over-expressed  in  cancer  (TF-BS)  were  cloned  in  tandem  upstream  of  a  minimal  lateADEp  promoter  (late 
adenovirus  promoter). 

(B)  Synthetic  tumor-specific  promoters  regulating  mKate2  exhibited  enhanced  tumor  specificity.  Representative 
native  tumor-specific  promoters,  SSXlp  and  H2Alp,  are  shown  for  comparison.  Parentheses  indicate  the 
transcription  factor  binding  sites  used  to  construct  each  synthetic  tumor- specific  promoter.  aHDF  (adult  human 
dermal  fibroblast)  and  HOV  (human  ovarian  epithelium)  are  normal  primary  cells.  OV8  (OVCAR8)  is  a  human 
ovarian  cancer  cell  line. 
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Figure  5.  Multi-output  AND  gate  architecture  and  performance. 

(A)  The  incorporation  of  a  synthetic  transcription  factor  (GAD)  as  the  output  of  Module  1  means  that  it  is  only 
expressed  at  high  levels  when  both  PI  and  P2  are  high  (state  [1,1]).  GAD  enables  facile  tuning  of  downstream 
gene  expression  (mKate2  •  ••  Output  N). 

(B)  mKate2  expression  can  be  modulated  by  modifying  the  number  of  GAD  binding  sites  in  the  synthetic 
GALp  promoter  expressing  mKate2  (G5,  G8,  G14  for  5,  8,  14  GAD  binding  sites  in  GALp)  and  the  presence  of 
miRNA  binding  sites  at  the  mKate2  3’UTR  (G5-Pe,  G8-Pe,  G14-Pe  indicate  transcripts  with  a  perfect  match 
miRl  binding  site).  Note  the  circuit  diagram  depicts  the  version  with  a  perfect  miRl  binding  site  in  the  mKate2 
transcript.  State  [0,0]  indicates  cells  containing  the  negative  control  modules.  Module  Icon  and  Module  Icon 
(see  text  for  details).  State  [0,1]  indicates  cells  containing  Module  Icon  and  Module  2.  State  [1,0]  indicates  cells 
containing  Module  1  and  Module  Icon.  State  [1,1]  indicates  cells  containing  Module  1  and  Module  2.  All  cell 
states  tested  also  contained  their  respective  Module  3. 

(C)  The  G8-Pe  AND  gate  architecture  triggered  high  mKate2  expression  in  human  ovarian  cancer  cells  (OV8) 
but  not  in  normal  cells  in  state  [1,1]. 

(D)  Our  circuit  design  prevents  potential  off-target  effects  in  primary  human  T  cells.  The  S(E2Fl)p  promoter 
was  active  in  T  cells  but  S(cMYC)p  was  inactive  in  T  cells.  The  G8-Pe  circuit  triggered  minimal  mKate2  output 
in  T  cells,  thus  minimizing  the  potential  off-target  effects  that  would  have  been  observed  if  the  S(E2Fl)p  was 
used  to  directly  drive  mKate2  expression. 

(E)  Our  circuit  design  can  be  readily  adapted  to  distinguish  a  different  tumor  type  from  its  normal  counterparts. 
By  using  two  different  promoters,  the  G8-Pe  circuit  triggered  a  strong  output  in  the  breast  cancer  cell  line 
MDA-MB-453  but  minimal  output  in  the  non-tumorigenic  breast  epithelial  cell  line  MCF-IOA. 

Error  bars  represent  S.E.M.,  N  =  3  biological  replicates  for  all  experiments  except  N  =  6  biological  replicates 
for  OVCAR8  group  in  Figure  5C. 


See  also  Figure  S4A. 


Figure  6 


25000-. 
20000- 
^  15000- 
^  10000- 
5000- 
0” 


STE  staining 


-T" 


STE: 


<5“'  <?-■  K*^ 

&  sN' 


B 


STE-mediated  lysis 


IFN-y  secretion 


IL12  secretion 


80000 
60000- 
S  -lOOOO- 

<>l 

^  20000- 


"W 


»Cs 


G 


anti-PD1  secretion 


0 


Figure  6.  Synthetic  circuits  triggered  human  ovarian  cancer-specific  expression  of  STE,  T-cell  mediated  lysis 
and  IFN-y  secretion  on  tumor  cells.  aHDF  (adult  human  dermal  fibroblast)  and  HOV  (human  ovarian 
epithelium)  are  normal  primary  cells. 

(A)-(C)  The  circuit  in  Figure  5A  was  engineered  to  display  STE  on  cell  surfaces  as  the  output.  STE  was 
encoded  in  Module  3.  State  ‘[1,1]  cont.’  indicates  cells  containing  Module  1  and  Module  2  and  a  negative 
control  Module  Scon  that  expresses  the  non-specific  transcription  factor  rtTA3  as  the  output.  State  [1,1] 
indicates  cells  containing  Module  1,  Module  2,  and  Module  3.  State  [1,0]  indicates  cells  containing  Module  1,  a 
negative  control  Module  2con  in  which  S(E2Fl)p  expresses  ECFP  without  miRl-BS(B),  and  Module  3.  Circuit- 
transduced  tumor  cells  but  not  normal  cells  expressed  high  levels  of  STE  in  state  [1,1],  which  triggered  robust 
T-cell  killing  and  IFN-y  secretion.  Student’s  t-test  was  performed  to  compare  the  [1,1]  states  between  normal 
cells  and  OV8  cells,  and  the  [1,1]  state  versus  the  [1,0]  state  in  OV8  cells. 

(D)-(G)  Synthetic  circuits  triggered  human  ovarian  cancer-specific  expression  of  combinatorial 
immunomodulators.  The  circuit  in  Figure  5A  was  extended  to  express  multiple  immunomodulatory  outputs 
(surface-displayed  STE  and  secreted  CCL21,  IL12,  and  anti-PDl  Ab)  specifically  in  OV8  cells  but  not  normal 
primary  cells.  State  ‘[1,1]  cont.’  indicates  cells  containing  Module  1  and  Module  2  and  a  negative  control 
Module  Scon  that  expresses  the  non-specific  transcription  factor  rtTA3  as  the  output.  State  [1,1]  indicates  cells 
containing  Module  1,  Module  2,  and  Module  3.  Each  immunomodulatory  gene  was  expressed  from  its  own  G8p 
promoter  encoded  on  a  lentivirus,  except  for  the  anti-PDl  Ab,  which  was  split  into  two  lenti viruses  encoding 
light  chain  (LC)  and  heavy  chain  (HC)  respectively.  Cell  lines  were  co-infected  with  these  lentiviral  constructs. 
Student’s  t-test  was  performed  to  compare  output  levels  in  the  [1,1]  states  between  normal  cells  and  OV8  cells. 

Error  bars  represent  S.E.M.,  N  =  3  biological  replicates  (*  p  <  0.05;  **  p<0.005). 

See  also  Figure  S4. 
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Figure  7.  Combination  immunotherapies  triggered  by  cancer- specific  circuits  significantly  reduced  tumor 
burden  and  increased  survival  in  an  intraperitoneally  disseminated  human  ovarian  cancer  model  in  mice. 

(A)  OV8  cells  were  transduced  with  constructs  encoding  a  control  output  rtTA3  (Control),  STE-only  (STE),  or 
combination  immunomodulators  (SCIP).  All  groups  were  implanted  into  NSG  mice  and  injected  periodically 
with  human  T  cells.  Student’s  t-test  was  performed  to  compare  tumor  burden  between  groups  at  day  37.  Error 
bars  represent  S.E.M.,  N  =  5  biological  replicates. 

(B)  Kaplan-Meier  survival  curves  of  various  groups  (N  =  5  biological  replicates).  A  log-rank  (Mantel-Cox)  test 
was  performed  to  compare  survival  between  groups. 

(C)  Robust  therapeutic  efficacy  was  achieved  even  when  only  15%  or  30%  of  tumor  cells  were  transduced  with 
the  SCIP  circuit  (the  remaining  implanted  cancer  cells  were  wild-type  OV8).  Therapeutic  efficacy  was 
dependent  on  the  presence  of  T  cells.  Student’s  t-test  was  performed  to  compare  tumor  burden  between  groups 
at  day  41.  Error  bars  represent  S.E.M.,  N  =  4  biological  replicates. 

(D)  Kaplan-Meier  survival  curves  of  various  groups  with  different  percentages  of  the  overall  tumor  cells 
transduced  with  the  SCIP  circuit  (N  =  4  biological  replicates).  A  log-rank  (Mantel-Cox)  test  was  performed  to 
compare  survival  between  groups. 

(E)  Lentiviral  delivery  of  the  SCIP-expressing  circuit  significantly  reduced  ovarian  cancer  burden.  NSG  mice 
were  injected  with  OVCAR8  cells  into  the  peritoneal  space  at  day  0.  Lentiviruses  encoding  circuit  Module  1 
and  Module  2  and  a  control  output  rtTA3  (Control),  or  combination  immunomodulators  (SCIP)  were  then 
injected  i.p.  on  day  7.  These  mice  were  also  injected  i.p.  periodically  (day  9,  16,  and  23)  with  human  T  cells. 
Student’s  t-test  was  performed  to  compare  tumor  burden  between  groups  at  day  39.  Error  bars  represent  S.E.M., 
N  =  5  biological  replicates. 

(F)  Kaplan-Meier  survival  curves  of  various  groups  (N  =  5  biological  replicates).  A  log-rank  (Mantel-Cox)  test 
was  performed  to  compare  survival  between  groups. 

(*p<0.05;  **p<0.005). 

See  also  Figure  S5-7. 
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Figure  SI.  Optimization  of  the  miRNA  binding-site  sequences  in  the  sponge  in  Module  2,  related  to  Figure  3. 

(A)  The  activity  of  the  promoters  SSXlp,  H2Alp,  and  CMVp  in  HEK-293T  cells. 

(B) -(C)  Various  miRNA  binding  sequences  (miRl-BSs)  used  in  the  sponge  affect  sponging  activity,  assayed  by 
mKate2  fold-activation  by  each  sponge  (state  [1,1]  versus  state  [1,0]  in  Figure  2B).  Nucleotides  highlighted  in 
dark  blue  are  mismatched  nucleotides.  Nucleotides  highlighted  in  red  are  insertions.  Light  blue  capital-delta  (A) 
indicates  deleted  nucleotides.  SSXlp  was  used  to  drive  Module  1  expression  and  H2Alp  was  used  to  drive 
Module  2  expression.  Module  1  encoded  a  transcript  with  miRl  as  an  intron  within  the  mKate2  gene,  where 
miRl  targeted  3  perfect  match  miRl-BS  encoded  downstream  of  mKate2.  Module  2  encoded  the  miRl  sponges, 
described  above,  downstream  of  an  ECFP  gene. 

(D)  MEI  of  mKate2  and  ECEP  in  binding  site  optimization  experiments.  Error  bars  represent  S.E.M.,  N  =  3 
biological  replicates.  MFI  =  Mean  Fluorescence  Intensity. 
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Figure  S2.  Optimization  of  the  sponge  architecture  in  Module  2,  related  to  Figure  3. 


(A)  The  number  of  miRNA  binding  sites,  the  sponge  expression  level,  and  the  architecture  of  the  sponge  affect 
sponging  activity.  Fold  activation  was  assayed  by  mKate2  fold-activation  by  each  sponge  (state  [1,1]  versus 
state  [1,0]  in  Figure  2B).  SSXlp  was  used  to  drive  Module  1  expression  and  H2Alp  or  CMVp  was  used  to  drive 
Module  2  expression,  as  indicated.  Module  1  encoded  a  transcript  with  miRl  as  an  intron  within  the  mKate2 
gene,  where  miRl  targeted  3  perfect  match  miRl-BS  encoded  downstream  of  mKate2.  Module  2  encoded  the 
miRl  sponges.  5B  and  lOB  refers  to  5  or  10  bulged  miRl  binding  sites  in  the  sponge  based  on  the  Bv2  design 
from  Figure  SIB.  ECl  and  EC2  correspond  to  exon  1  and  exon  2  of  ECFP  and  constructs  in  which  lOB  was 
encoded  within  the  ECFP  transcript  with  different  architectures  (Cons  and  ZRANBl).  12TuD  (Tough  Decoy) 
architecture  refers  to  24  miRl-BSs  encoded  within  12  loop  of  an  RNA  hairpin  (2  miRl-BSs  per  hairpin  loop). 
5B/100  bp  refers  to  5  bulged  miRl-BSs  separated  by  lOObp  linkers  rather  than  5bp  in  all  other  sponges.  Error 
bars  represent  S.E.M.,  N  =  3  biological  replicates. 

(B)  MFI  of  mKate2  and  ECFP  in  sponge  architecture  optimization  experiments. 
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Figure  S3.  pri-miRNA  backbones  used  to  produce  miRl  from  Module  1  affect  circuit  performance,  related  to 
Figure  3. 

(A)  Module  1  variants  Mvl-Mv3  expressed  a  transcript  in  which  the  different  miRl  pri-miRNA  architectures 
shown  were  encoded  as  an  intron  within  the  mKate2  gene,  where  miRl  then  targeted  3  perfect  match  miRl-BS 
encoded  downstream  of  mKate2.  The  colored  circles  denote  the  location  of  the  insertion  sequences  for  creating 
Mv2  and  Mv3  (See  Table  S2  for  detailed  sequences  and  functional  pri-miRNA  modifications). 

(B-E)  Different  pri-miRNA  backbones  modulate  the  production  of  miRNA  and  the  ON:OFF  ratio  of  the  circuit. 
(B-C)  We  tested  the  three  different  Module  1  variants  together  with  Module  2  variations  composed  of  a  control 
sponge  with  no  miRl  binding  sites  (‘Cont.’)  and  a  sponge  with  5  perfect  match  miRl  binding  sites  downstream 
to  an  ECFP  gene  (‘5Pe’).  Module  1  was  expressed  by  SSXlp  and  Module  2  was  expressed  by  H2Alp.  Mvl- 
Mv3  exhibited  different  miRl  production  efficiencies,  as  measured  by  ECFP  signal  reduction  with  the  5Pe 
Module  2  compared  to  ECFP  signal  in  the  control  Module  2;  Mv3  had  the  most  efficient  miRl  production. 

(D-E)  mKate2  expression  levels  with  control  sponge  or  the  optimized  sponge  (‘Sponge  7’,  Sponge  7  from 
Figure  S2A).  Mv3  exhibited  a  low  mKate2  level  in  the  OFF  state  and  the  highest  fold- activation  by  the  sponge 
in  the  ON  state. 

Error  bars  represent  S.E.M.,  N  =  3  biological  replicates. 
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Figure  S4.  Multi-output  circuit  specifically  triggered  T  cells  to  kill  ovarian  tumor  cells  and  secrete  IFN-y, 
related  to  Figure  5-6. 

(A)  FACS  histograms  showing  the  activity  of  single  promoters  and  the  G8-Pe  circuit  when  transduced  into  T 
cells  via  lentivirus.  The  percentages  on  the  histograms  denote  the  percent  of  cells  with  higher  mKate2 
expression  level  than  the  value  of  the  dotted  line.  A  constitutive  human  Ubiquitin  C  promoter  expressing 
mKate2  was  used  for  the  “infection  control”  group. 

(B)  S(cMYC)p,  S(E2Fl)p,  and  the  G8-Pe  circuit  exhibit  high  activity  in  OV8  cells  but  not  in  various  non- 
tumorigenic  cell  lines.  OV8  is  an  ovarian  cancer  cell  line.  HOF  (human  ovarian  fibroblasts),  HOMEC  (human 
ovarian  microvasculature  epithelial  cells),  CCD-841 -CoN  (non-tumorigenic  colonic  epithelial  cells),  and  MCF- 
lOA  and  MCE-12A  (non-tumorigenic  mammary  epithelial  cells),  are  non-tumorigenic  cell  lines.  iPSCs  = 
inducible  pluripotent  stem  cells. 

(C-E)  T-cell  mediated  killing  of  G8-Pe  circuit-transduced  tumor  cells  was  significantly  greater  than  of  normal 
cells.  The  percentage  of  cell  death  mediated  by  T  cells  from  donor  #1  (D),  donor  #2  (E)  and  the  calculated  mean 

(C)  are  shown.  *  indicates  condition  not  tested.  Note  that  total  cell  death  and  STE-mediated  cell  death  of  the 
[1,0]  state  in  aHDF  and  HOV  cells  were  not  measured  because  this  state  simulates  a  condition  in  which  only 
one  of  the  input  promoters  is  active  and  is  expected  to  result  in  lower  killing  than  state  [1,1],  in  which  both 
input  promotes  are  active.  Total  cell  death  and  STE-mediated  cell  death  correspond  to  uncorrected  cell  death 
and  background-killing  corrected  cell  death,  respectively.  See  Methods  for  more  details. 

(F)  G8-Pe  circuit-transduced  tumor  cells  trigger  robust  IFN-y  secretion  by  T  cells.  Data  shown  represent  EFN-y 
levels  secreted  by  T  cells  from  donor  #1  and  donor  #2.  Error  bars  represent  S.E.M.,  N  =  3  biological  replicates. 

(G-I)  aHDF,  HOV,  and  OV8  constitutively  displaying  STE  can  trigger  robust  STE-mediated  killing  and  IFN-y 
secretion  by  T  cells.  A  constitutive  human  Ubiquitin  C  promoter  (hUbCp)  was  used  to  drive  STE  expression  on 
aHDF,  HOV,  and  OV8.  All  three  cell  lines  displayed  STE  and  triggered  strong  STE-mediated  cell  killing  and 
EFN-y  secretion  by  T  cells.  Error  bars  represent  S.E.M.,  N  =  3  biological  replicates. 

(J)  Doxycycline  (Dox)-inducible  expression  of  STE  triggered  significant  T-cell-mediated  killing  in  vivo. 
Inducible  display  of  STE  with  doxycycline  (Dox)  on  the  surface  of  OV8  cells  triggered  human  T-cell-mediated 
killing  in  vivo  in  a  mouse  model  of  disseminated  intraperitoneal  ovarian  cancer.  Anti-tumor  activity  was 
abrogated  when  STBs  were  not  displayed  (T  cells  only)  or  T  cells  were  not  injected  (Dox  only).  Student’s  t-test 
was  used  to  compare  tumor  burdens  between  groups  at  day  45.  Error  bars  represent  S.E.M.,  N  =  5  biological 
replicates 

(K)  Tumor  burden  of  individual  mice  (black)  and  mean  (red)  of  the  Dox  +  T  cell  group  in  Figure  S4J. 

(L)  OV8  cells  engineered  with  a  Dox-inducible  STE  construct  expressed  high  levels  of  STE  on  their  cell 
surfaces  after  induction.  Histograms  of  STE  staining  on  OV8  cells  are  shown.  The  gray  and  orange  histograms 
represent  the  staining  of  uninduced  OV8  and  Dox-induced  OV8  cells,  respectively.  Data  shown  are 
representative  histograms,  which  were  repeated  in  3  independent  experiments. 


(*p<0.05;  **p<0.005;  ***p  <0.001). 
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Figure  S5.  Combination  immunotherapies  triggered  by  the  circuits  significantly  reduced  tumor  burden  in  an 
intraperitoneally  disseminated  ovarian  cancer  model  within  NSG  mice  periodically  injected  with  human  T  cells, 
related  to  Figure  7. 

(A)  Tumor  burden  of  individual  mice  (black)  and  mean  (red)  of  each  tested  group  (N  =  5  biological  replicates). 
X-axes  denote  days  post  tumor  injection,  Y-axes  denote  tumor  burden  (photons/sec). 

(B)  In  vivo  tumor  bioluminescence  images  of  each  mouse  from  day  5  (D5)  to  37  (D37).  This  data  corresponds 
to  the  experiment  shown  in  Figure  7A-B. 
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Figure  S6.  Robust  therapeutic  efficacy  was  achieved  even  only  15%  of  tumor  cells  were  transduced  with  the 
SCIP-expressing  circuit  and  then  implanted  into  the  peritoneal  cavity  of  mice,  related  to  Figure  7. 

(A)  Tumor  burden  of  individual  mice  (black)  and  mean  (red)  of  each  tested  group  (N  =  4  biological  replicates). 
X-axes  denote  days  post  tumor  injection,  Y-axes  denote  tumor  burden  (photons/sec). 

(B)  In  vivo  tumor  bioluminescence  images  of  each  mouse  from  day  5  (D5)  to  41  (D41).  This  data  corresponds 
to  the  experiment  shown  in  Figure  7C-D. 
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Figure  S7.  Lentiviral  delivery  of  SCIP  eireuit  mediated  robust  therapeutie  effieaey  in  an  intraperitoneally 
disseminated  human  ovarian  eaneer  model  in  miee,  related  to  Figure  7. 

(A)  Tumor  burden  of  individual  miee  (blaek)  and  mean  (red)  of  eaeh  tested  group  (N  =  5  biologieal  replieates). 
X-axes  denote  days  post  tumor  injeetion,  Y-axes  denote  tumor  burden  (photons/see). 

(B)  In  vivo  tumor  biolumineseenee  images  of  eaeh  mouse  from  day  13  (D 13)  to  39  (D39).  This  data 
eorresponds  to  the  experiment  shown  in  Figure  7E-F. 
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Table  SI.  Comparison  of  native  and  synthetic  promoter  sequences.  We  analyzed  the  sequences  of 
representative  native  (H2Alp  and  SSXlp)  and  synthetic  S(cMYC)p  and  S(E2Fl)p  cancer-specific  promoters 
with  Matinspector  software  to  predict  transcription  factor  binding  sites  (TF-BSs).  Total  matches:  total  number 
of  TF-BSs  identified;  TF  families  (annotated  as  ‘Matrix  families’  by  Matinspector):  a  group  of  different  TFs 
that  share  very  similar  TF-BSs  (i.e..  Position  Weighted  Matrix)  and  biological  function;  Average  matches  per 
family:  the  mean  number  of  TF-BSs  per  matrix  family,  calculated  by  the  total  number  of  identified  TF-BSs 
divided  by  the  number  of  different  TF  families;  Max  single-family  matches:  the  maximum  number  of  TF-BSs 
that  belonged  to  a  single  TF  family;  Max  family  %  of  total:  the  percentage  of  ‘Max  single-family  matches’  out 
of  the  total  number  of  TF-BSs  identified;  Matches  per  lOObp:  the  mean  number  of  TF-BSs  per  lOObp  in  the 
promoter,  used  to  normalize  TF-BSs  by  promoter  length.  Bottom  -  graphical  depiction  of  TF-BS  families  on 
promoter  sequences.  Tags  indicate  the  location  and  orientation  of  each  TF-BS  (above  and  below  the  line  for 
forward  and  reverse  strands,  respectively).  Each  different  color  denotes  a  different  matrix  family. 
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ITAAGTGTGCT  CGCTTC  GGCAGCACATATAC  TATGTT  GAATGAG  GCTTCAGTACTTT  AC  AGAA 
;TAAGTGTGCT  CGCTTC  GGCAGCACATATAC  TATGTT  GAATGAG  GCTTCAGTACTTT  AC  AGAA 


Mvl  TCGTTGCCTGCACATCTTGGAAACACTTGCTGGGATTACTTCTTCAGG- 
Mv'2  TCGTT  GCCTGC ACATCTT  GGAAAC ACTTGCT  GGGATTACTTCTT  CAGG - 


— TTAACCCAACAGAAGGCTCGAG - TGCTG 

— TTAACCCAACAGAAGGCTCGAQAAGGTATATTGCTG 

MV'S  TCGTTGCCTGCACATCTTGGAAACACTTGCTGGGATTACTTCTTCAGGGACTTCTTAACCCAACAGAAGGCTCGAGAAGGTATATTGCTG 
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Mvl  GA - GAATTCAAGGGGCTACTTTAGGAGCAATTATCTTGTTTACTAAAACTGAATACCTTGCTATCTCTTTGATACATT 

MV'2  GACrrCAAGGGGCTAGAATTCAAGGGGCTACTTTAGGAGCAATTATCTTGTTTACTAAAACTGAATACCTTGCTATCTCTTTGATACATT 
Mv3  GACrrCAAGGGGCTAGAATTCAAGGGGCTACTTTAGGAGCAATTATCTTGTTTACTAAAACTGAATACCTTGCTATCTCTTTGATACATT 
Insert  1  EcoRI 
+CNNC  motif 


TTTACAAAGCTGAATTAAAATGGTATAAATTAAATCACTTTTTTCAATTGTACTAACTTCGAG^flmmgmmmmgjmmmmm^CACAGl 

tttacaaagctgaattaaaatggtataaattaaatcacttttttcaattgtactaacttcgagZS5S55E5S55EScacag| 

TTTACAAAGCTGAATTAAAATGGTATAAATTAAATCACTTTTTTCAATTGTACTAACTTCGAGtiaUitMiiUiUiUiUlUitMilUiUlSCACAGl 
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Table  S2.  Design  and  optimization  of  the  miRl  expression  cassette  in  Module  1,  related  to  Figure  3. 

(A)  Structure  of  pri-miRNA  in  Module  1.  Structural  features  of  pri-miRNAs  important  for  efficient  pri-miRNA 
processing  and  frequently  observed  features  in  mammalian  pri-miRNAs  are  indicated:  basal  UG  motif,  35  bp 
double  stranded  stem  structure,  CNNC  motif  in  the  3’  flanking  region,  apical  UGU/GUG  motif  in  terminal  loop 
and  mismatch  at  position  8. 

(B)  Sequence  alignment  of  pri-miRNA  designs  (Mvl,  Mv2  and  Mv3). 
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Supplementary  Information 

Sponge  optimization  in  Module  2 

We  first  optimized  the  miRl  binding-sites  (miRl-BSs)  in  Module  2  for  maximal 
sponging  efficiency.  We  designed  14  variants  of  miRl-BSs  with  different  mismatches,  insertions, 
and  deletions  compared  to  perfect-match  miRl-BSs  (Figure  SIB-D).  We  measured  the  sponging 
efficiency  of  each  variant  by  comparing  mKate2  levels  generated  by  Module  1  in  the  presence  of 
the  Module  2  sponge  encoding  5  identical  tandem  miRl-BSs  in  the  3’  UTR  of  ECFP  versus  no 
miRl-BSs  (hereafter  reported  as  “fold  ON-OFF  induction”).  As  an  indication  of  miRl  binding  to 
the  sponge,  we  also  expected  to  observe  a  decrease  in  ECFP  expression  in  constructs  where  the  3’ 
UTR  of  ECFP  contained  miRl-BSs  versus  no  miRl-BSs.  We  used  two  highly  active  native 
promoters  in  HEK-293T  cells,  SSXlp  (Cure  et  al.,  1997)  and  H2Alp  (Rogakou  et  ah,  1998),  to 
express  Module  1  and  Module  2,  respectively.  In  HEK-293T  cells,  H2Alp,  which  expresses  the 
sponge  in  Module  2,  is  >3-fold  stronger  than  SSXlp,  which  expresses  Module  1  (Figure  SI  A). 
Thus,  the  Module  2  sponge  transcript  concentrations  were  expected  to  be  in  excess  compared  to 
the  Module  1  transcript,  which  ensured  enhanced  sponging  efficiency.  When  we  tested  a  sponge 
transcript  with  perfect  match  to  miRl  generated  by  Module  1,  we  found  a  80%  decrease  in  the 
levels  of  ECFP,  even  though  only  a  33%  decrease  was  expected  due  to  the  1:3  ratio  of 
miRTsponge  transcript  concentrations.  This  phenomenon  may  reflect  recycling  of  miRNA 
molecules  following  target  transcript  inhibition,  which  is  especially  efficient  for  perfect-match 
binding  sites  (Baccarini  et  al.,  201 1)  (Figure  SID). 

We  calculated  the  binding  efficiency  of  miRl  to  its  corresponding  binding  sites  in  the 
sponge  (A  AG)  with  a  previously  described  model  (Kertesz  et  al.,  2007),  but  we  found  that  A  AG 
did  not  predict  sponging  efficiency  (Figure  SIB).  Specifically,  calculations  for  microRNA 
binding  efficiencies  to  target  sites  were  performed  at: 
http://genie.weizmann.ac.il/pubs/mir07/mir07  prediction.html.  For  example,  binding  sites  with 
perfect  complementary  to  miRl  (Pe)  had  the  strongest  binding  efficiency  (AAG=  -26),  but 
limited  mKate2  activation  (Figure  SIC,  1.7-fold  ON-OFF  induction),  which  could  be  a  result  of 
efficient  sponge  transcript  degradation  due  to  perfect  miRNA  complementarity  (Ebert  and  Sharp, 
2010).  A  AAG  =  -15  was  calculated  for  several  binding  sites  that  contained  various  mismatches 
to  miRl  that  resulted  in  bulged  complementarity  to  miRl  (Bv2,  Bv4,  Bv6  and  Bvl5).  Yet,  the 


31  measured  ON-OFF  induction  for  these  binding  sites  ranged  from  1.25-fold  (Bvl5)  to  2.62-fold 

32  (Bv2),  which  covered  almost  the  entire  dynamic  range  of  mKate2  activation.  miRl  binding  site 

33  version  Bv2  contained  mismatches  to  nucleotides  9-12  of  miRl  and  exhibited  the  best  sponging 

34  efficiency  (2.6-fold  mKate2  ON-OFF  induction),  likely  because  mismatches  at  these  positions 

35  protect  the  sponge  from  Ago2-mediated  degradation  while  allowing  efficient  miRl  binding  by 

36  preserving  seed  sequence  complementarity  (Ebert  and  Sharp,  2010).  We  therefore  used  Bv2  for 

37  the  sponge  in  further  optimization  of  the  circuit.  Hereafter,  we  use  the  abbreviation  ‘#B’  to  refer 

38  to  the  number  of  tandem  Bv2  binding  sites  separated  by  5  bp  spacers  (for  example,  ‘5B’  refers  to 

39  5  Bv2  binding  sites).  Finally,  we  also  observed  that  ECFP  levels  could  not  predict  mKate2 

40  activation.  These  results  highlight  the  complexity  of  miRNA-mediated  regulation,  in  which  the 

41  complementarity  of  the  miRNA  to  its  cognate  binding  sites  determines  not  only  the  inhibitory 

42  pathway  to  which  the  target  mRNA  is  designated  (mRNA  cleavage,  destabilization  or  translation 

43  inhibition),  but  also  the  rate  of  unloading  of  miRNAs  from  Ago2,  which  is  enhanced  by  high 

44  complementarity  of  the  miRNA  to  its  target  (De  et  al.,  2013).  While  bulged  BSs  generally  have 

45  better  efficiency  in  shunting  miRNAs  compared  to  perfect  match  BSs,  we  believe  that  deriving 

46  optimized  rules  on  how  to  introduce  effective  bulges  will  require  additional  studies. 

47  We  then  tested  how  the  number,  spacing,  location  and,  architecture  of  miRl  Bv2  binding 

48  sites  in  the  sponge  transcript,  together  with  sponge  expression  levels,  affect  sponging  efficiency. 

49  We  observed  2.9  and  3.6-fold  mKate2  ON-OFF  induction  when  5B  versus  lOB  binding  sites 

50  were  encoded  in  the  3’  UTR  of  the  Module  2  sponge  transcript  expressed  by  H2Alp, 

51  respectively  (Figure  S2A,  Sponge  1:  H2Alp-ECFP-5B  versus  Sponge  4:  H2Alp-ECFP-10B). 

52  Switching  the  promoter  in  the  sponge  encoding  the  5B  binding  sites  from  H2Alp  to  CMVp,  a 

53  4.3-fold  stronger  promoter  compared  to  H2Alp  (Figure  SI  A),  resulted  in  a  6-fold  ON-OFF 

54  induction  in  mKate2  levels  (Figure  S2A,  Sponge  3:  CMVp-ECFP-5B).  Therefore,  increasing  the 

55  sponge  transcript  concentration  efficiently  increases  miRl  shunting,  suggesting  that  the  promoter 

56  selected  to  express  Module  2  should  be  stronger  than  the  one  selected  to  express  Module  1. 

57  Switching  the  location  of  the  lOB  from  the  3’  UTR  to  the  5’  UTR  of  ECFP  in  Module  2  {Sponge 

58  4:  H2Alp-ECFP-10B  versus  Sponge  5:  H2Alp-10B-ECFP)  reduced  the  sponging  efficacy 

59  from  3.6-fold  mKate2  ON-OFF  induction  to  2.4-fold  (Figure  S2A),  which  suggests  that  the 

60  miRl-BS(B)  were  more  accessible  when  encoded  in  the  3’  UTR  of  the  mRNA  rather  than  in  the 

61  5’  UTR.  Interestingly,  increasing  the  spacing  between  miRl-BSs  located  in  the  3’  UTR  of  the 


62  sponge  transcript  from  5  bp  to  100  bp  {Sponge  2:  H2Alp-ECFP-5B/100bp)  only  mildly  reduced 

63  the  sponging  efficiency  from  2.9  to  2.3-fold  mKate2  ON-OFF  induction.  Additionally,  encoding 

64  24  miRl-BSs  in  the  loop  of  an  RNA  hairpin  with  previously  described  tough-decoy  (TuD) 

65  architecture  (2  miRl-BSs  per  loop,  12  loops  overall.  Sponge  9:  H2Alp-ECFP-12TuD  (Ebert 

66  and  Sharp,  2010;  Haraguchi  et  ah,  2009)  provided  only  3.8-fold  mKate2  ON-OFF  induction 

67  (Figure  S2A). 

68  Finally,  to  increase  the  molar  ratio  of  sponges  produced  from  a  given  promoter,  we 

69  designed  a  Module  2  variant  where  lOB  were  encoded  as  an  intron  within  an  ECFP  coding 

70  sequence  in  addition  to  5B  encoded  in  the  3’UTR  of  the  transcript  (Figure  S2A).  In  this 

71  architecture,  the  donor  and  acceptor  consensus  sequences  that  flank  the  10  intronic  miRl-BS(B) 

72  were  derived  from  a  previously  described  synthetic  intron  {Sponge  6:  H2Alp-ECl-10Bcons- 

73  EC2-5B)  (Greber  et  ah,  2008).  Based  on  ECFP  expression,  this  sponge  was  correctly  spliced, 

74  thus  producing  ECFP  mRNA  that  contained  5B  and  an  intron  with  lOB.  This  sponge 

75  demonstrated  ~4-fold  mKate2  ON-OFF  induction,  slightly  better  than  Sponge  4:  H2Alp-ECFP- 

76  lOB. 

11  We  designed  an  additional  architecture  in  which  the  intronic  consensus  sequences  that 

78  flank  the  lOB  were  derived  from  the  ZRANBl  gene  (Memczak  et  ah,  2013)  {Sponge  7:  H2Alp- 

79  EC1-10Bzranbi-EC2-5B).  This  design  was  supposed  to  generate  a  stable  circular  lOB  exon  and 

80  EC1/EC2-5B  introns,  thus  disrupting  ECFP  expression.  This  architecture  demonstrated  enhanced 

81  sponging  efficacy  with  ~6-fold  mKate2  activation  (Figure  S2A-B)  and  did  not  produce  ECFP 

82  signal.  To  test  whether  that  the  lack  of  ECFP  signal  was  caused  by  the  expected  improper 

83  splicing,  resulting  in  an  impaired  ECFP  transcript,  or  because  the  transcript  was  not  spliced  at  all, 

84  we  generated  an  additional  sponge  in  which  the  ZRANBl  splicing  sequences  were  deleted 

85  {Sponge  8:  H2Alp-ECl-10B-EC2-5B).  Sponge  8  demonstrated  almost  identical  sponging 

86  efficacy  as  Sponge  7  (Figure  S2A),  which  supports  the  hypothesis  that  Sponge  7  did  not  splice. 

87  Since  Sponge  7  and  Sponge  8,  which  did  not  splice,  demonstrated  higher  sponging  efficiency 

88  than  Sponge  6,  which  did  splice,  we  concluded  that  intronically  encoded  miRNA  sponges  did  not 

89  significantly  contribute  to  sponging  activity.  Thus,  Sponge  7:  H2A1p-EC1-10Bzranbi-EC2-5B 

90  architecture  was  chosen  as  the  final  ‘optimized  sponge’  design  for  our  circuit. 


91 


92  miRNA  backbone  optimization  in  Module  1 

93  We  hypothesized  that  our  circuit  performance  could  be  optimized  by  enhancing  miRl 

94  production  from  Module  7.  We  therefore  modified  the  pri-miRNA  stems  for  efficient  pri-miRNA 

95  processing  (Figure  S3A  &  Table  S2)  (Auyeung  et  ah,  2013;  Fellmann  et  ah,  2013;  Suzuki  et  ah, 

96  2015).  Previous  studies  have  described  multiple  features  of  pri-miRNA  stems  that  promote 

97  Drosha-mediated  pri-miRNA  processing,  including  a  UG  motif  at  the  base  of  the  pri-miRNA 

98  hairpin,  a  narrow  range  of  tolerable  pri-miRNA  stem  lengths,  a  mismatched  motif  in  the  basal 

99  stem  region,  and  a  CNNC  motif  16-18bp  downstream  to  the  Drosha  processing  site  that  is 

100  required  for  efficient  SRp20-dependent  processing  of  the  pri-miRNAs  (Auyeung  et  ah,  2013; 

101  Fellmann  et  ah,  2013).  Our  miRl  design  Mvl  was  based  on  a  previously  described  intronically 

102  encoded  synthetic  miRNA  (Greber  et  ah,  2008;  Nissim  et  ah,  2014;  Xie  et  ah,  2011)  derived 

103  from  a  miR-30  backbone  (Fellmann  et  ah,  2013).  The  £'co7?/ restriction  site  at  position  231-236 

104  was  generated  by  mutating  the  original  sequence  of  endogenous  miR-30  from  GACTTC  to 

105  GAATTC  (where  position  1  is  the  first  intronic  nucleotide  following  the  splicing  of  Module  1). 

106  This  point  mutation  disrupted  an  essential  CNNC  motif  in  the  miRNA  backbone  and  resulted  in 

107  reduced  pri-miRNA  processing  efficiency.  The  Xhol  restriction  site  at  position  128-133  was 

108  generated  by  mutating  the  original  sequence  of  the  endogenous  miR-30  from  CTAAAG  to 

109  CTCGAG,  which  could  further  impair  the  pri-miRNA  secondary  structure. 

110  We  generated  our  Mv2  design  by  adding  a  CNNC  motif  (CTTCAAGGGGCTA) 

111  downstream  to  base  pair  231  of  Mvl  {Insert  1)  and  an  additional  insert  (AAGGTATAT) 

112  downstream  to  base  pair  133  of  Mvl  {Insert  2)  to  improve  the  secondary  pri-miRNA  structure 

1 13  that  was  altered  by  the  formation  of  the  Xhol  restriction  site  (Fellmann  et  ah,  2013)  (Figure  S3A 

114  &  Table  S2).  To  generate  our  Mv3  design,  we  added  an  additional  insert  (GACTTC) 

115  downstream  to  base  pair  111  of  Mvl  {Insert  3)  to  improve  the  secondary  pri-miRNA  structure 

116  that  was  altered  by  the  formation  of  the  EcoRI  restriction  site  (Fellmann  et  ah,  2013)  (Figure 

117  S3  A  &  Table  S2).  This  resulted  in  a  pri-miRNA  backbone  similar  to  the  optimized  miR-E 

118  backbone,  which  is  superior  to  the  miR-30  backbone  (Fellmann  et  ah,  2013). 

119  We  compared  miRl  production  from  each  backbone  by  characterizing  their  efficiency  at 

120  inhibiting  mKate2  levels  in  the  Module  I  transcript,  in  which  they  were  encoded,  and  ECFP 

121  levels  in  a  Module  2  mRNA  containing  5  perfect  match  miRl-BSs  (5Pe)  (Figure  S3B-E).  As 


122  reference,  we  used  a  control  construct  with  no  miRl-BSs  (Figure  S3B).  Pri-miRNA  Mv3 

123  mediated  a  93%  decrease  in  ECFP  expression  in  the  Module  2  construct  with  the  5Pe  sponge 

124  compared  to  the  control  construct,  which  was  superior  to  Mvl  (80%  inhibition)  and  Mv2  (71% 

125  inhibition),  demonstrating  that  pri-miRNA  design  can  tune  the  efficiency  of  miRl  production. 

126  This  data  suggests  that  that  miRl  production  was  more  potent  in  Mv3  compared  to  Mv2  and 

127  Mvl. 

128  We  also  validated  the  effects  of  the  miRl  backbone  on  the  efficiency  of  mKate2 

129  production  by  measuring  mKate2  activation  levels  in  the  presence  of  the  optimized  sponge  in 

130  Module  2  (Figure  S3D-E).  The  absolute  mKate2  signal  in  the  OFF  and  ON  states  (e.g.  mKate2 

131  levels  measured  with  control  versus  optimal  sponge)  were  both  affected  by  the  pri-miRNA 

132  design:  Mvl  exhibited  an  OFF-ON  shift  from  -1400  to  -8400  Mean  Fluorescence  Intensity 

133  (MFI)  (5.8-fold  activation  ratio),  Mv2  from  -170  to  -500  MFI  (3-fold  activation)  and  Mv3  from 

134  -380  to  -2600  MFI  (6.8-foId  activation).  Based  on  the  ECFP  repression  data  from  comparing  the 

135  5Pe  sponge  to  the  control  construct,  we  expected  Mv3  to  be  more  efficient  at  inhibiting  mKate2 

136  than  Mvl  and  Mv2  in  the  presence  of  the  5Pe  sponge.  However,  Mv2  exhibited  potent  inhibition 

137  of  mKate2  levels  in  both  the  OFF  and  ON  states.  This  suggests  that  our  modifications  to  the  pri- 

138  miRNA  design  may  not  only  affect  the  efficiency  of  miRl  production,  but  also  other  regulatory 

139  mechanisms.  We  hypothesize  that  these  pri-miRNA  modifications  could  modify  the  efficiency  of 

140  existing  splicing  sites  or  may  have  added  alternative  splicing  sites  in  the  miRNA  intron  that 

141  could  impair  mKate2  production  by  alternative  splicing.  This  hypothesis  is  supported  by  the  fact 

142  that  the  pri-miRNAs  in  designs  Mv2  and  Mv3  include  a  potential  splicing  donor  sequence 

143  AAGGT  introduced  by  Insert  2,  and  that  a  potential  donor  site  CAGGT  in  Mvl  and  Mv2  was 

144  modified  to  CAGGG  by  Insert  3  in  Mv3  (Table  S2). 

145  As  noted  in  the  main  text,  we  selected  the  Mv3  architecture  for  further  experimentation 

146  since  it  had  low  mKate2  level  in  OFF  state  and  exhibited  the  highest  fold-activation  by  the 

147  sponge  in  ON  state. 

148  Finally,  the  miRNA  duplex  sequence  can  be  modified  to  promote  asymmetric  guide 

149  strand  selectivity  that  favors  miRl  production.  The  miRl  pri-miRNA  is  processed  into  a  miRNA 

150  duplex  capable  producing  both  miRl  encoded  by  the  duplex  3p  arm  and  an  alternative  miRNA 

151  encoded  by  the  duplex  5p  arm  (Table  S2).  The  ratio  of  3p  (miRl):5p  production  is  determined  by 


152  the  superposed  patterns  of  the  5'-end  nucleotide  identity  and  the  thermodynamic  stability  of  two 

153  miRNA  duplex  termini  (Suzuki  et  al.,  2015).  Quantitative  thermodynamic  and  sequence  analysis 

154  of  the  miRl  duplex  has  shown  that  the  3p  arm  has  5'-uridine  and  a  thermodynamically  unstable 

155  5'-end,  whereas  5p  arm  has  5'-guanine  and  a  thermodynamically  stable  5’-end  (Suzuki  et  al., 

156  2015).  This  indicates  that  the  miRl  duplex  is  optimized  to  produce  mostly  the  3p  arm  (e.g. 

157  miRl);  thus,  it  was  therefore  not  modified.  Nevertheless,  we  note  that  asymmetric  guide  strand 

158  selectivity  is  an  important  consideration  in  designing  synthetic  miRNAs. 

159 

160  Bioinformatics  framework  for  synthetic  promoters  design 

161  We  identified  genes  that  were  differentially  expressed  in  human  ovarian  cancer  cells 

162  compared  to  healthy  tissues  from  the  publicly  available  The  Cancer  Genome  Atlas  (TCGA) 

163  database  (Bell  et  al.,  201 1)  using  R  (version  3.0.2).  We  followed  standard  procedures  to  calculate 

164  differentially  expressed  genes  between  cancer  cases  compared  to  controls  using  t-test.  The  p- 

165  values  were  corrected  for  multiple  testing  using  the  Bonferroni  correction.  We  then  selected 

166  genes  that  were  previously  identified  as  transcription  factors  (TFs)  from  the  literature 

167  (Vaquerizas  et  al.,  2009;  Wingender  et  al.,  2013).  From  these  genes,  we  selected  the  ones  that 

168  were  hypomethylated  in  cancer.  For  that  purpose,  we  used  the  publicly  available  database  of 

169  DNA  Methylation  and  gene  expression  in  Human  Cancer  MethHC  (Huang  et  al.,  2015).  This 

170  database  provides  the  most  hyper-  and  hypomethylated  genes  for  different  cancers.  We  selected 

171  TF  genes  with  the  most  hypomethylated  promoter  regions  using  default  settings.  TF  expression 

172  in  OVCAR8  cells  was  confirmed  using  the  CCLE  database  (Barretina  et  al.,  2012).  We  then 

173  determined  the  binding  motifs  for  ovarian  cancer-enriched  TFs  using  GREAT,  MEME,  JASPAR 

174  and  MOTIFMAP  (Bailey  et  al.,  2009;  Daily  et  al.,  2011;  McLean  et  al.,  2010;  Sandelin  et  al., 

175  2004)  and  assembled  them  into  our  synthetic  promoters  (Figure  4A-B). 

176 

177  Analysis  of  Natural  and  Synthetic  Promoters 

178  We  analyzed  promoter  sequences  using  Matinspector  software  to  predict  transcription 

179  factor  binding  sites  (TF-BSs)  (Cartharius  et  al.,  2005).  We  ran  our  analysis  on  two  natural 

180  promoters  (SSXlp,  H2Alp)  and  two  synthetic  promoters  (S(E2Fl)p,  S(cMYC)p)  using  core 

181  similarity  0.9  and  matrix  similarity  optimized  +  0.05.  We  examined  individual  TF-BSs  as  well  as 


182  TF-BS  families  (defined  as  TFs  with  highly  similar  TF-BS  -  for  example,  the  binding  sites  for 

183  MYC,  USF  and  Max  were  grouped  as  the  family  E-BOX).  We  found  that  the  synthetic 

184  promoters  had  more  predicted  TF-BSs  per  lOObp.  For  example,  S(cMYC)p  contained  35.8  TF- 

185  BS  matches  per  100  bp  while  H2Alp  and  SSXlp  contained  ~7.9  TF-BS  matches  per  lOObp.  The 

186  synthetic  promoters  had  matches  to  less  TF-BS  families  (~6-7  different  TF-BS  families) 

187  compared  with  the  natural  promoters  (-40-41  different  TF-BS  families).  The  synthetic  promoters 

188  had  more  TF-BS  matches  for  each  of  the  different  TF-BS  families  compared  with  the  natural 

189  promoters.  Furthermore,  we  identified  the  TF-BS  family  with  the  highest  number  of  matches  for 

190  each  promoter,  and  found  that  the  top  TF-BS  family  for  the  synthetic  promoters  had  more  TF-BS 

191  matches  than  the  natural  promoters.  For  example,  E2F  TF-BS  family  was  found  to  be  66%  of  the 

192  total  predicted  TF-BSs  in  S(E2Fl)p  and  the  E-BOX  TF-BS  family  (corresponding  to  MYC,  USF, 

193  MAX)  were  found  to  be  29%  of  the  total  predicted  TF-BSs  in  S(cMYC)p.  In  contrast,  the  ZF5F 

194  (ZFP  161)  TF-BS  family  was  found  to  be  only  11%  of  the  total  predicted  TF-BSs  in  H2Alp 

195  while  the  MYBL  (MYB)  TF-BS  family  was  found  to  be  only  6%  of  the  total  predicted  TF-BSs 

196  in  SSXlp. 

197 

198  Multi-output  AND  Gate  optimization  for  Combinatorial  Immunomodulation 

199  We  first  characterized  the  output  levels  of  four  distinct  circuit  states  in  OVCAR8  cells 

200  (Figure  5B).  To  do  so,  we  encoded  each  module  on  separate  lentiviruses  and  infected  OVCAR8 

201  cells  with  the  corresponding  Module  3  architecture  and  either:  (A)  a  negative  control  Module  1 

202  that  expresses  the  non-specific  transcription  factor  rtTA3  {Module  Icon)  and  a  negative  control 

203  Module  2  in  which  S(E2Fl)p  expresses  ECFP  without  miRl-BS(B)  {Module  2con)  to  model  the 

204  scenario  when  neither  cancer-specific  promoters  are  active  [0,0];  (B)  Module  Icon  and  Module  2 

205  to  model  the  scenario  when  only  S(E2Fl)p  is  active  [0,1];  (C)  Module  1  and  Module  2con  to 

206  model  the  scenario  when  only  S(cMYC)p  is  active  [1,0];  (D)  Module  1  and  Module  2  to  model 

207  the  scenario  when  both  promoters  are  active  [1,1]. 

208  The  circuits  expressed  minimal  mKate2  fluorescence  in  states  [0,0]  and  [0,1].  In  state 

209  [1,0],  there  was  a  low  background  mKate2  level  that  increased  from  5  to  8  GAL4-BSs  in  the 

210  synthetic  GALp  and  was  attenuated  with  the  presence  of  miRl-BS(Pe)  in  Module  3  (Figure  5B). 


211  Only  when  both  input  constracts  were  present  (state  [1,1])  did  we  observe  a  significantly 

212  increased  mKate2  output.  Furthermore,  the  mKate2  output  level  could  be  enhanced  by  increasing 

213  the  number  of  GAL4-BSs  in  the  Module  3  synthetic  promoter  from  5  to  8  GAL4-BSs.  Although 

214  adding  perfect  match  miRl  binding  sites  in  the  3’  UTR  of  the  mKate2  transcript  decreased  the 

215  basal  and  maximal  output  levels,  it  also  improved  fold-activation  (G5-Pe,  G8-Pe,  and  G14-Pe 

216  versus  G5,  G8,  and  G14,  respectively).  The  G8-Pe  circuit  exhibited  negligible  mKate2  levels  in 

217  states  [0,0]  and  [0,1],  and  high  8-fold  mKate2  activation  in  state  [1,1]  over  the  minimal  mKate2 

218  expression  in  state  [1,0],  and  was  thus  chosen  for  further  use. 


